CERN-PH-TH/2010-089 



The warped, resolved, deformed conifold 

gets flavoured 

Jerome Gaillard°[^ Dario Martelli^ Carlos Nunez'j^ loannis Papadimitriou^ 

a: Department of Physics, Swansea University 
Singleton Park, Swansea SA2 8PP, United Kingdom 

b: Department of Mathematics, King's College, London 
Strand, London WC2R 2LS, United Kingdom 

c: Department of Physics, CERN Theory Unit 
1211 Geneva 23, Switzerland 



Abstract 
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conifold with additional flavour fivebrane sources. The resulting solution is a generalisation 
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D3 brane sources. We propose that this solution may be interpreted in terms of a combined 
effect of Higgsing and cascade of Seiberg dualities in the dual field theory. 
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1 Introduction 

The systematic study of supersymmetric geometries of String/M theory has provided new 
valuable tools for addressing a number of problems in the context of the gauge/string duality 
as well as in string compactifications. Powerful techniques for studying supersymmetric 
backgrounds of various supergravities are based on G-structures and generalised geometry. In 
this paper we will consider a class of supersymmetric backgrounds of Type IIB supergravity 
characterised by an SU{3) structure [2], and we will discuss a simple solution generating 
method applicable to these geometries. In particular, we will show that starting from a 
solution to the "torsional superstring" equations of [3], a more general interpolating solution 
may be generated, that includes the simple class of warped Calabi-Yau solutions P] in a limit. 
In most cases, this procedure is equivalent to the chain of U-dualities discussed in [5j, as was 
also showed in [B]. However, exploiting the relation to generalised calibrations, this method 
can be applied as well to geometries which include the back-reaction of supersymmetric 
sources. In fact, we will apply the procedure to a supersymmetric solution describing Nc D5 
branes wrapped on the S"^ inside the resolved conifold, plus Nf D5 branes sharing the M}'^ 
Minkowski directions and infinitely extended along a transverse cylinder [7j. We will study 
the case in which the Nf sources are smeared over their transverse compact directions and 
we can have Nf/Nc ~ C>{1). The motivations and consequences of this smearing have been 
discussed in detail in [S]. 

The "seed" solution on which the solution generating technique mentioned will be applied 
was discussed in [9rl The large radius (UV) behaviour of this solution is such that the dilaton 



asymptotes to a constant In section [3.2.2 we will study the small radius (IR) asymptotics 
suitable for the purposes of this work, extending the analysis in [HI [ID]- The solution we will 
discuss is singular at the origin of the radial direction and we may view this singularity as 
one that can be resolved, in a fashion similar to the Klebanov-Tseytlin solution [11], where 
away from the singularity the solution captures the correct Physics. We will not resolve this 
singularity in the present paper, but we anticipate that one way to do this is to consider a 

-"^More precisely, the solutions in [9] had different IR asympotics to the ones we will present in this paper. 
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profile for tlie smeared flavour branes that vanishes smoothly at the origin [12]. The solution 
depends on two integers Nc and Nf, and we will fix the boundary conditions in the IR by 
requiring that setting Nj = 0, the smooth solution discussed in |5j is recovered. 

The field theory proposed to be dual to the solutions in [3, [HI HU] is a deformation of A/" = 1 
SQCD by coupling the quarks to the infinite tower of massive states (Kaluza-Klein modes) 
whose presence is due to the twisting of the original six- dimensional field theory. See [13] 
for details about the twisting and the spectrum. 

After applying the solution generating technique to the "seed" solution, we will find a 
new background with non-zero Ramond-Ramond (RR) and Neveu-Schwarz (NS) fields. In 
addition to the integers Nc, Nf, the solution depends on three continuous parameters. One 
of them is the string coupling at infinity. Then there is a parameter we denote c, which is 
related to the size of the 5*^ as measured from infinity [3] , hence to the amount of resolution 
in the geometry. In the unfiavoured case this parameter was related to the VEV of baryonic 
operators in the baryonic branch of the Klebanov-Strassler theory [Tl]. Finally, we have a 
parameter introduced by the transformatiorj^ The solution may be viewed as a "flavoured" 
version of the resolved deformed conifold solution [1], thus justifying our title. Sending to 
zero the resolution parameter (c — oo), and Nf, we will obtain a solution closely related to 
the Klebanov-Strassler solution. 

We will then discuss the dual field theory interpretation of our new solutions. We begin 
with the simpler solution mentioned above: in this case we will argue that the presence of 
smeared source D3 branes may be interpreted as a change of ranks of the gauge theory due 
to Higgsing, that in our case will take place at all energy scales. This is inspired by the ideas 
discussed in [T2|. More precisely, we will propose a Klebanov-Strassler type quiver theory 
SU {Nc + n + Uf) X SU {n + Uf) in the mesonic branch, where the change of ranks is due to 
a combinatiorl^ of Seiberg dualities (running of n) and of Higgsing (running of Uf). 

The field theory interpretation of the general case when we restore a non-zero value of Nf 
is more complicated. We will make some general comments and then we will analyse two 
kinds of field theory scenarios. In both cases we will consider a quiver gauge theory with 
two nodes, each with different ranks of the gauge groups, and we will assume that there are 
bi-fundamental fields interacting with a quartic superpotential, precisely as in the Klebanov- 



Strassler theory. In section 4.2.1 we will analyse a quiver in which we have added explicit 
flavours, i.e. fields transforming in fundamental representations of the gauge groups. In this 
case the global flavour symmetry is manifest along the flow from the UV to the IR. However, 



^As we will explain later, this parameter corresponds to the boost parameter in [S]. 

■^We are indebted to Ofer Aharony for crucial comments that helped us sharpening this proposal. 
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in this picture a precise understanding of the interplay of Seiberg duahties and Higgsing 



appears problematic, as we will explain. In section 4.2.2 we will discuss the possibility that 
the flavour symmetry may develop dynamically after Higgsing and the cascade take place. 
As we will see, in this case we can obtain a rather detailed matching between field theory 
and gravity computations. 

The rest of the paper is organized as follows: in section [2] we will discuss a solution gener- 
ating method applicable to supersymmetric Type IIB geometries characterised by an SU (3) 
structure. We will also explain how to incorporate supersymmetric sources (flavour branes). 
In section [3] we will give details of the new solutions. First we will extend previous studies 
on fivebrane solutions with asymptotically constant dilaton, and then we will construct ex- 
plicitly the new flavoured resolved deformed conifold solution. In section |4] we will discuss 
possible dual field theory interpretations of these backgrounds. The two scenarios alluded to 
above will be discussed in detail. Section [5] contains a discussion of a Z2 symmetry, and its 
breaking, in gravity and in the field theory. Finally, we summarise our findings in section [6j 
The appendices include the Type IIB equations of motion and various technical calculations 
that we will quote in the paper. 

2 Generating solutions from SU(3) structures 

We start this section by presenting the BPS equations of an SU{3) structure background, 
derived from the general set-up in [2]. While the authors of [2] work with pure spinors, 
in a particular case, their results can be formulated in terms of the two differential forms 
characterising the SU{3) structure. First, let us state precisely the ansatz and conventions. 
We work with Type IIB supergravity in Einstein frame and consider a ten-dimensional 
space as a warped product of a four-dimensional Minkowski space and a six-dimensional 
space equipped with an SU (3) structure, using the metric 

ds^ = e^'^[dxl, + dsl]. (2.1) 

We also have several fluxes: the RR forms F(3), F(5) and the NS three- form H. These 
have components only in the internal six-dimensional space, except for F(5) that is self-dual. 
Generically, we can write 

F(5)=e^^+*(l + *io)vol(4)A/, 
H = dB, 
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where / is a one-form. The general Type IIB supersymmetry conditions for these geometries 
were derived in [2], [16] as equations for the two pure spinors (multi-forms) \E'i, \E'2 and read 

/ <f,\ ,-„A+5<I'/4 

(d-iJA)[e2^-*/2v[/2] = 0. 

(2.3) 

We then specialise these to the case of SU{3) structure. This means that the two pure 
spinors take the form 

^1 = -—e^+'^/^ ( 1 - ze2^+*/V _ Ae^A+f J ^ j\ = q (2.4) 

8 y 2 J 8 

Here J is the (would-be) Kahler two-form and Q is the holomorphic three-form. Together, 
these define an SU{3) structure on the six- dimensional geometry. In addition, we have a 
function ( arising as a phase in the pure spinor \E'i. For this reason, the SU{3) structure is 
referred to as interpolating 

Equating terms involving forms of the same degree we obtain the BPS equations of the 
system written as 

d {e^^+^/'n) = 

d(e«^JA J) =0 
d(e2^-*/2cosC) =0 

-e-^^-*d(e^^sinC)=/ (2.5) 
-e* cos C *6 i^{3) - e2^+3*/2 ^.^ ^.^ ^) ;\ J = e-2A-*/2d (e4A+* J) 

- sin Ce* *6 i^(3) + cos Ce'^+^^^/M (e"* sin C) A J = 

-^d(e-*sinC) AJAJ = *6i^(i) 

Manipulating these equations a little more one can show that 

if = d(tanCe2^+*/V) ^ 5 = tanCe^^+^/V , (2.6) 

thus the (non-closed part of the) B field is determined by the SU{3) structure. Notice 
that the first equation in (2.5) implies that the geometry is complex, in the usual sense, as 
opposed to the general case discussed in [2], [16]. This gives a useful characterisation of the 
geometries we are interested in. These equations were also derived in which discussed 
first the results presented in this section. 

^See [T71 [TH] for earlier work on interpolating geometries. 
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In the rest of this paper, we will impose that = 0. The last equation in (2.5) then 
implies that d (e~* sin (^) = and the system simplifies further, reducing to 



d (e^^+^/^fi) =0, d (e2^-*/2 cos C) = 0, d {e'^J A j) = 0, 



— e 



cosC *6 F(3) = e-2^-*/M (e^^+* J) , H = - sinCe* H ^(3), 
sinC) = /• 



(2.7) 



It is instructive to speciahse the system (2.7) to the case C = 0: 

d (e6^+*/2fi) =0, d (e2^-*/2) =0, d {e'^J A J) = 



-e* *6 ^(3) 



g-2A-*/2^ (e^^+* J) 

/ = 0, = 0, = 0. 



The only non-zero fiux is then F(3) and the BPS system describes a configuration of D5 
branes. These are simply the S-dual version of the "torsional superstring" equations of 
[31 Hn] and they were written in this form in [201 El] • A notable solution to these equations 
was discussed in 



We will now show that from a solution of the system (2.8) one can generate a solution 



of the more complicated system (2.7) for a non- vanishing (. This is then a simple solution 
generating technique. We will sometimes refer to this procedure as rotatior^ Precisely, we 
have that if one defines 



Q 
J 

^(3) 
"(5) 



$(0) 
1 



2A(0) 



cosC 
cosC 



^1 ^$/2 

cosC 



J(0) 



(2.9) 



-(1 + *io)dx° A dx^ A dx'^ A dx^ A d 



sinC ^4^(0) 



cos^ C 



where the quantities with a ^ ' obey the equations in (2.8), then the (new) quantities on the 
left-hand side obey the equations in (2.7). ki is here an integration constant. We require 



^There are different motivations for this: firstly, it is a rotation in the space of KiUing spinors, parame- 
terised by C- Secondly, in the particular case discussed in [5] this corresponds to an actual rotation of NS5 
branes in a T-dual Type IIA brane picture. 
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that the condition = is preserved, which imphes that d (e *sin(^) = 0. We can then 
solve this equation obtaining 



sin^ = 



(2.10) 



where K2 is another integration constant. This formula requires the dilaton to be bounded 
from above at any position in space. To summarise, let us write the background after the 
"rotation" in terms of the initial one (in Einstein frame). 



B 



-<I>/2 



**6d(e2*jW) 

J(0) 



(0)2 
6 



(2.11) 



-^2(1 + *io)vol(4) A dh ^ 



where 



2 (e 



-2'S> 



4) 



(2.12) 



Any solution of the system (2.5), supplemented by the Bianchi identities for the fluxes, is 



a solution of the equations of motion of Type IIB supergravity [21 |T6]. One can then show 



that imposing the Bianchi identities for the simplified (seed) system (2.8), implies also the 
Bianchi identities, and hence the full equations of motion, of the more complicated system. 



Thus, starting from a solution of the system (2.8) (with a bounded dilaton), one can generate 



a solution of the system (2.7) using the formulas in (2.9). This result was discussed also in 



2.1 Adding D5 brane sources 

In this subsection we will show how the generating technique discussed above may be ap- 
plied also to supersymmetric solutions for a combined system of supergravity plus smeared 
sources. The key observations are the following. Firstly, when the sources are smeared in 
a supersymmetric way, the D-brane action can be written in terms of generalised calibra- 
tions and their net effect is captured by simple modifications of the Bianchi identities for 
the fluxes. In particular, the non-closed part of the fluxes is identified with the so-called 
smearing form - see [23]. The supersymmetry equations, however, do not change in form. 
Then by using the results of [21] , the computation of the previous subsection can be applied 
to the case with sources. The interest of including the back-reaction of such explicit branes 
is that these may be interpreted as flavours in the context of the gauge/gravity duality. 
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We consider the case of D5 brane sources. One then needs to study the combined action of 
Type IIB supergravity with the DBI and WZ terms for the source branes. Using the SU{3) 
structure cahbration conditions, the combined action can be written as [23] 

S = Sub - j (e^^+*/Vol(4) A (cosCe^^J + sinCe-*/^^) - C(6) + C(4) A A S(4) (2.13) 

where Sub is the action of Type IIB supergravity. Here 2(4) is the smearing form, charac- 
terising the distribution of sources. It is proportional to the number of flavour branes Nj 
and has no components along the Minkowski directions. C(4) and C(6) are defined via 



^(5) 



dC(4) + 5 AF(3), 



e* *ioF(3) =dC(6)+5AF(5). 



(2.14) 



Using the results of [21] , we know that the addition of sources, even when smeared, does not 



modify the form of the BPS system (2.5) but only the Bianchi identities. These now read 

(2.15) 



dF(3) = S(4), 
dF(5) = HA F(3) 



5 AS 



(4). 



As shown in 



if one imposes the Bianchi identities, every solution of the BPS system 



is a solution of the equations of motion coming from (2.13). As previously described, one 



can generate a solution to the equations of motion of the action (2.13) from the case C = 0- 



Setting = in (2.13) and using the fact that 5 = we find 



D5 sources 



j (e6^+*/2^ol(4)AJ-C(6)) AH(4). 



(2.16) 



This is the action for supersymmetric D5 brane sources in a background characterised by 



the equations (2.8) 



2.2 The limit of D3 brane sources 

The limit C — > 7r/2, in which the supergravity background goes over to the warped Calabi-Yau 
geometry, is slightly more subtle [5j. Here we will determine how the action for the source 
branes behaves in this limit, and we will find that indeed the limiting action corresponds to 
smeared source D3 branes, with a particular smearing form arising in the limit. Considering 



the action (2.13) as generated from the C = case (2.16), we can work out the dependence 



on C, of every quantity from eqs. (2.7) and (2.9 ). Recalling that the RR potentials are defined 
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using (2.14), we have 

.2A 



— 1 2A(0) '^l 'I>/2 

e = -e = -t ' 

cos C cos C 

$ 1 . ^ 

e = — sm 



S = ^^sin3/2Ce'^'"'jW 
^6) = ^Vv^e6^*"Vol(4) A J(°) 

^1 V ^^2 



(2.17) 



sinC 4A(o) , 

W4) = ^0l{4) 

cos^ C 

Looking at the dependence on cos C, the action for the sources can be written as 

5'sourccs = ^j^([cos2CsinCe^^vol(4) A5 + cos2CC(4) AS] +0(cos'C)) A 2(4). (2.18) 

In this formula, the quantity in square brackets does not scale with cos^, implying that this 
is approaching a finite non-zero value when goes to 7r/2, and additional terms in 0(cos^ Q 
go to zero in the limit. However, there is an overall factor cos~^ C,. Therefore, if we want the 
action to be finite in the limit cos^ — ?■ 0, then we need to scale S(4) accordingly. However, 
S(4) = A^/W(4) where a;(4) does not depend on C,. We then conclude that we need to impose 
the following condition: 

— — )■ constant when C, ^ —. (2-19) 



cos2 C 2 
In this case the limit of the D5 brane source action is 



5'd3 sources = - / ( 6^"^'° Vol(4) + C(4) ) A 2(6) (2.20) 



1(4) ^ L.(4) j l\ ^(6) 

where the tilded quantities correspond to the limit of the untilded ones. We have defined 
H(6) as 

B A 2(4) S(6) when C ^ 2 ' '^^'^^^ 
In the limiting case we can then identify the action as the smearing of supersymmetric D3 
branes with smearing form S(6). 

3 Adding D3 branes to the flavoured D5 brane solution 

Here we apply the procedure discussed in the previous section to a solution representing D5 
branes wrapped on the S"^ of the resolved conifold, with the addition of explicit smeared D5 
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brane sources. The resulting Type IIB solution with D3 brane charge and B field will be a 
"fiavoured" version of the warped resolved deformed conifold solution originally derived in 
[T]. We will postpone a field theory interpretation of this new background until section |4| 

3.1 D5 branes on the resolved conifold with flavour D5 branes 

The setup corresponding to D5 branes wrapped on the S*^ of the resolved conifold, with 
addition of smeared D5 sources was described in [71ini[ID|. The metric in the Einstein frame 
takes the form 

ds^ = e^^ [d43 + e*-^^dsg , (3.1) 

where here the internal metric dsg does not change under the "rotation" procedure. As we 
have seen in the previous section, the un-rotated metric is obtained by setting e^^ = e^^*"' = 



KiC*/^ in (3.1). The solution can then be completely describecjj in terms of the sechsbeins 
e", a = 1, . . . , 6, parameterising the internal metric dsg, i.e. dsg = 6abe°'e^- In the notation 
of [7] the sechsbeins are 



= ^e^ui + aui) , = ^e^ (w2 - a^a) , = ^e'=(c:;3 + ^3), (3.2) 
where the one-forms Ui, a)j, i = 1, 2, 3, are defined as 

ui = d9, = cos '?/'d6' + sin T/^ sin ^d(^, 

L02 = sm9dip, a;2 = — sin'?/'d6' + COS'?/' sin6'd(^, (3.3) 
W3 = cos 9dip, 0J3 = dip + cos 9dip. 

Moreover, the RR three-form is given by 

F(3) = -2Ar,e-23"VAe2Ae3 + ^(a2-2a6+l-^)e-2''-VAe^Ae=^ 
+Nc{b - a)e-^-''-\e^ A e'^ + A e^) A 

+ ^b'e-9-'i-''eP A (-e^ A + e'^ A e^). (3.4) 

As was shown in [9], there is a set of variables that decouples the BPS equations described 
in the previous section for this ansatz and leads to a single second order ordinary differen- 
tial equation, whose solution completely determines the supergravity background. See also 
Appendix |D} This set of variables is introduced defining 

e''=l{pS^), e^^ = P cosh r-g, e^>^ = AY, 

^ Psinhr u cr (3 5) 



P cosh T-Q' Nc' 

^In Appendix [d| we have written this solution in the variables used in [S]. 
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Solving the resulting set of decoupled BPS equations one then finds 

1 



sinhr 



sinh(2(p - po)) 



Q = {Qo + 



2Nr - N. 



coshr 



2Nr - N. 



(2p coshr — 1) 



a = tanh r ( Q 



2Nr-N 



(3.6) 



,4($-$o) 



cosh^(2po) 



Y=-{P' + Nf), 

while the only remaining unknown function, P{p), is determined by the equation 

P" + iP' + N,) (^^^4^ + ^^^^ - 4cosh.l ^ 0. (3.7) 



P-Q 



P + Q 



Here po, Qo and are constants of integration and we set Qo = —Nc + Nf/2. The SU{3) 
structure for this class of backgrounds is specified by a (would-be) Kahler form J and a holo- 
morphic three-form Q, which may be written explicitly as follows (cf. for the unfiavored 



case 



j(o) ^ A + A (-cos/ie"^ + sin/ie^) + A (-sin/ie'^ - cos/ie^) , (3.8; 
^(0) ^ (e'' -Fie^) A [e^ + z (-cos/^e'^ -Fsinpe^)] A [e^ +2 (-sinpe'^ - cos/ie^)], 

where the angle < p < 7r/2 corresponds to a rotation in the e*^ — plane and is given by 

P — coshrQ 



cos /i 



P cosh T — Q 



(3.9) 



The SU{3) structure for the transformed solution is now obtained simply from (3.8) via the 



rescalings (2.9). 



Charge quantisation 



Given a solution of (3.7) one immediately obtains the full string background via the above 



relations. In particular, the new background is obtained as in (2.11) and it depends on the 



parameters Nc and Nf, which in the case C = 0; can be interpreted respectively as the number 
of colour and fiavour D5 branes. However, this interpretation should be reconsidered for the 
transformed backgrounds. Notice that in the presence of sources one should be careful with 
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Fs\Nf=o = N,T,, (3.10) 



the D5 charge quantisation condition for the original background. In particular, since F3 
is not closed, its integral over the three-cycle at infinity will depend on the representative 
submanifold, hence it cannot be quantised. We therefore define the number of colour D5 
branes by integrating over the three-cycle F3 evaluated at Nf = 0. Since the latter is closed, 
this definition makes sense and we have 

1 

where is any representative of the unique three-cycle at infinity. For the transformed 
background this should be modified to 

/ F3U,=o = N,T, = ^T, e N. (3.11) 
2^10 J S3 «i 

We also redefine the number of fiavours using 

dFg = ^ sin ^ sin OdO A d(/? A d^ A d^ = ^ sin ^ sin OdO A d(/? A d^ A d^. (3.12) 
4 Aki 

Noticing that P, Q and J^^'^ are homogeneous of A*",, and Nf of degree one, it follows that 
the rotated background, in string frame, becomes 

dsi, = h-^/Mxls + e^''h^/^K^dsl{N„Nf) 
F3 = e-2* *6d(e2*j(°)(iV„iV;)), 
S = «:2e2*j(0)(iV„iV;), 

i^5 = -/«2(l + *io)d/i"^Avol(4). (3.13) 



Using the expression (2.12) for h, it is clear that the constant ki can be absorbed into the 
rescaled charges N^ and Nf. Namely, defining 

h = - kI (3.14) 

and absorbing ki by a trivial rescaling of the worldvolume coordinates — )■ we have 

dsl,. = h-'/^dxl, + e'^h'/^dsliN,, Nf) 
F3 = e-2* *6d(e2*j(°)(iV„iV^)), 
B = K2e'''J^'\N,,Nf), 

i^5 = -/«2(l + *io)d/i"^Avol(4). (3.15) 

It follows that the effect of the rotation described in the previous section is simply the 
introduction of the parameter K2, with < K2 < max{e^*}. For K2 = we recover the 
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original background. To make contact with the discussion in [5], we may parameterise K2 as 



fi;2 = e-*°°tanh/3, (3.16) 

where $00 is the asymptotic value of the dilaton. In [5] this transformation was derived 
as a simple chain of U-dualities, and the constant /3 arose as a boost parameter in eleven 
dimensions. However, the derivation presented here (see also |6]) may be readily applied to 
cases with sources. 

3.2 The flavoured resolved deformed conifold 

We will now present a deformation of the solution of Butti et al pQ, describing the baryonic 
branch of the Klebanov-Strassler theory [H], induced by the back-reaction of source D5 
branes. We will start by first reviewing some of the material in [5]. 

3.2.1 Review of the unflavoured solution 

Before we present the flavoured solution, let us recall the unflavoured solution [H [TJ 13 E]. 



This solution is obtained by setting Nf = in (3.15 ), Qo = —Nc, Po = and picking a speciflc 



solution of the differential equation (3.7) for P. The solution is only known numerically, but 
one can easily determine its IR and UV asymptotic forms, which are specifled in terms of 
two arbitrary constants, hi and 

' h n -U 4/11 fl 47V2\ 3 16/11 f . 47V2 32N^\ 5 r)( n^^ n ^ 



15 lif J P ^ 525 \ 3/if 3hf J 

P=<( (3.17) 



In the full solution the two constants are related in a non-trivial way [S], given in (B.41) 
of Appendix [B} What is important for the present discussion is that hi (c) as a function of 
c takes values in [2,-|-cxd), while c G [0, +00), with hi{0) = 2. One can also construct the 
solution in an expansion for large c, as is discussed in detail in Appendix |Bj One then flnds. 



via (|3.6|), that the dilaton takes the form 

1 

c 



^ ^2$^ ( ^ _ _hKs{p) + 0{l/c') ] , (3.18) 



''These constants are related to the parameters 7, too and C/ in as follows: 

D C 

Moreover, Nc here is M in |5J. The complete map to the variables used in [5 includes: = 2p, tmm = t, 

cmm = P/Nc, and /mm = 4:V/Nc- See also Appendices [b| and [p} 
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where 
and 

hKs = 2^1' Nl — (V coth(2p') - 1) (sinh(4pO - Ap'f" , (3.20) 

is the Klebanov-Strassler warp factor (cf. eq. (90) in [H]). 

Let us recall some limits of this two-parameter family of solutions discussed in [5] : 

• /3 ^ 

This is the original background before adding the D3 brane charge [7j, namely it 
describes wrapped D5 branes. The interpretation of the parameter c was discussed in 
Taking c — )■ corresponds to going to the near-brane limit, which is the solution 
discussed in [22]. In this decoupling limit the theory on the fivebranes was argued to 
flow to pure A/" = 1 in the IR [22]. In the opposite c — )■ oo limit, the metric approaches 
the deformed conifold with three-form flux. 

• /3 oo 



This limit ensures that the constant term in the warp factor in (3.14) is removed and 
the leading term in the UV is dominated by hxs- The expansion in large c does not 
terminate and c remains as the only non-trivial parameter of the solution in addition 
to $oo- This solution describes the baryonic branch of the Klebanov-Strassler theory 
and the parameter c is related to the baryonic branch VEV as ?7 oc c^^ [I]. 

/3 — )■ oo and c — )■ oo 



In this case the warp factor h in (2.11 ) is replaced by hxs and hence the background is 
the Klebanov-Strassler background [H]. In particular, the unwarped internal metric is 
the deformed conifold metric. The only free parameter is the asymptotic value of the 
dilaton $oo- The deformation parameter e of the deformed conifold may be reabsorbed 
by a rescaling of the metric. 

/3 — )■ oo and c — )■ 

This is the limit of large VEVs on the baryonic branch solution. In [S] it was shown that 
there is a large region where the dilaton is approximately constant and the solution 
is well approximated by the solution of [25]. It was then argued that in this limit 
one approaches the wrapped fivebrane theory, but with a B field on the two-sphere. 
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induced by the "rotation" procedure. Therefore in this case the theory is well described 
by fivebranes wrapped on a fuzzy two-sphere. 



3.2.2 The new flavoured solution 



Let us now present the new flavoured resolved deformed conifold solution. Note that the 
analysis so far is general enough to allow for smeared D5 sources and therefore we only 
need to find a new solution of the "master equation" (3.7) with A''^ 7^ 0, subject to the 



condition that it reduces to the solution of [T] in the limit Nf — )■ 0. We have found this new 
solution numerically, but again one can systematically determine the IR and UV asymptotics, 
discussed in detail in Appendix [B] 



P 



hp + ^ (^-p log p- log p) + ^ [ 



p log (-log p) 



logp 



+ C>(p3 log p), p^O, 



4p/3 I 9^/ I 1 



13\ 81^1 



-4p/3 



0{pe 



-8p/3^ 



P — 7- 

(3.21) 

These solutions are different from those discussed in |9] and in particular they reduce to 
their unfiavoured counterparts in (3.17) in the limit Nf — )■ 0. However, contrary to the 



unfiavoured solution, the flavoured solution is singular in the IR. Moreover, note that the 
flavours, i.e. terms proportional to Nf, dominate the IR, as well as the UV when /3 — cxd 



(see (3.24) below). Via (3.6) we see that this asymptotics imply that the dilaton goes to —00 
in the IR and not to a constant as in the unfiavoured case. In particular, for small values of 
the radial coordinate we hav43 



^2* 



N 



3e2*-c3/2 



1 



log (-logp) 
Slog p 



>/ (-logp)^/^ 

Thus the warp factor has the following expansion in the IR 

3/2 

(-logp)^/^ 



O 



logp 



(3.22) 



h 



2A^ 



/ 



1 + 



log (-logp) 



8 logp 



logp 



(3.23) 



The numerical solution interpolating between the asymptotic behaviours in (3.21 ) is plotted 



in Fig. [T| In Fig. |2] the numerical solution is explicitly compared with the asymptotic 



solutions (3.21) by zooming in the IR and UV regions. 



'To recover the limit Nf = at fixed c the expansions (3.22) and (3.23) are not usefuL This limit is 



completely smooth as can be seen from the expansions of P in (3.21 1 and the various plots of the numerical 
solutions. 
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Figure 1: Plot of the function P{p) and the dilat on for the numerical solution interpolating 
between the two asymptotic behaviours in (3.21). The plots correspond to the values N,. = 
10, Nf = 20, c = 30, and = 100. 
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Figure 2: In these plots we plot the same numerical solution as in Fig. [T} but we zoom in 
on the IR region (left) and on the UV region (right) and we compare the numerical solution 



(black) with the corresponding asymptotic solutions given in (3.21). These are plotted in 
red (IR solution) and in blue (UV solution). 



Let us now reconsider the various limits discussed in the previous section for the un- 
flavoured solution: 



/3 ^ 

Again, this is the original D5 brane background before adding the D3 branes. Further 
taking c — )■ is the near-brane (decoupling) limit. In this case, the resulting solution 



interpolates between the IR asymptotics given in (3.21) and the linear dilaton asymp- 



totics P ~ \1N^ - Nf\p + Po (with = ^)- The solution is plotted in Figs. |3|and|4| 
and it is the flavoured generalisation of the wrapped D5 solution of [22j . 

/3 oo 

Although in this limit we remove the constant term from the warp factor in the UV, 
now the leading form of the warp factor is not dominated by hxs any more, but by 



16 



the term introduced by the sources. In particular, using the expression (B.24) for the 
dilaton we have 

(o2/3j\T roo \ 
y dp' (sinh(4p') - 4p')"'^' + 0{l/c^) \ . (3.24) 

It follows that the UV asymptotic behaviour of the flavoured solution is different from 
the Klebanov-Strassler asymptotics. As we will see later, this leads to a different field 
theory picture. 

• c — )■ oo and /3 — t- cxd 

This limit cannot be taken naively in the flavoured case, due to the fact that the flavours 
dominate the UV after the leading constant term in the warp factor is removed. The 
reason why we cannot go to the Klebanov-Strassler limit, while keeping the flavour 
D5 branes, is that this is not a supersymmetric configuration |26]. Therefore we will 
consider the limit of c — )■ oo and Nj — )■ at fixed cNf. This limit will be the subject 
of the next subsection. In section |4| we will argue that the field theory interpretation 
of this solution is a modification (by Higgsing) of the Klebanov-Strassler cascade. 

• /3 — 7- oo and c — 

One can perform an analysis similar to that in [5] for the solution in this range of 
parameters and show that there is again a large region where the solution is well 
approximated by a resolved conifold metric, with addition of fluxes and sources. This 
suggests that there should exist an exact Type IIB solution analogous to that in [25|, 
modified by the presence of sources. It would be interesting to find this solution. 

Summary 

Let us summarise the effects of the addition of the flavour D5 branes to the unflavoured 
solution. The ansatz for the metric and fluxes is essentially unchanged and may be param- 
eterised completely in terms of the function P - see Appendix [D] for a concise presentation 
of the ansatz. The function P for various values of Nj is plotted on the left in Figure 5 
Notice that the leading UV behaviour of P is not affected by the flavours. However, the 
flavoured solution is actually singular at p = 0. After introducing the D3 branes, the six- 



dimensional metric dsg is unchanged, and is warped by the warp factor (3.14). This picks 
up the sub-leading behaviour of the dilaton in the UV and therefore it is sensitive to the 
Nf flavours - see the plots of h on the right of Figure 5 The divergence of /i at p = is 
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Figure 3: Plot of the function P{p) and t he di 



between the IR asymptotic behaviour in (3.21) and the hnear UV asymptotics P 

Pn 



aton for the numerical solution interpolating 



\2Nr 



Nf\p. The plots correspond to the values Nc = 20, Nf 



20, 



10 and hi = 25.93. 



Contrary to the Nj = case, we do not have an analytic expression for the value of hi 
leading to linear dilaton asymptotics in the UV. 



1.0 r 



P(P) 
250 - 



0.0 I — ^ 
0.000 



Figure 4: Here we plot the same numerical solution as in Fig. [3} but we zoom in on the IR 
region (left) and on the UV regio n and we compare the numerical solution (black) with the 
IR asymptotic solution given in (3.21) and the UV asymptotic solution P ~ \2Nc — Nf\p. 



These are plotted in red (IR solution) and in blue (UV solution) 



due to the singularity in the IR. The fall-off at infinity is noticeably slower with respect to 
the unflavoured case, and we expect that this will persist after resolving the IR singulatity. 
To understand the physical origin of the UV behaviour, we will next discuss the solution 
in a limit in which the six-dimensional (unwarped) metric becomes an ordinary deformed 
conifold. 



3.3 Adding smeared D3 branes to the Klebanov-Strassler theory 



We now discuss a solution obtained in the limit c — )■ oo (with tanh /3 = 1) of the flavoured 
solution. This limit can be obtained by inserting the warp factor (3.24) in ( |3.14 ) and sending 
c — i- oo. However, the fact that /i ~ 1/c and not /i ~ 1/c^, as is the case for the unflavoured 
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Figure 5: On the left: plots of P{p) for fixed values c = 30, Nc = 10 and different values 
of Nf (and hi). The continuous curve is Nj = 0. Superimposed on this are the curves for 
the following values: Nf = 5 (dotted green), Nf = 10 (dotted red), Nf = 20, (dotted blue), 
Nf = 40 (dotted black). On the right: different plots of h{p) for the same values of Nf. 

solution, does not allow to take this limit directly. To obtain a well-defined limit we set 



cNf = u, and keep u fixed in the limit c — ?■ oo. Using (B.24), we obtain 



h 



(^2^/=^// dp' (sinh(4p') - 4p')"'^' + hKs^ + (^(l/c')- (3.25) 



Inserting this in (3.15) and sending c — )■ oo we obtain an exact solution (notice that the 
expressions for B and below do not scale with the parameter c) 

d4. = e*- [h-'/'dxl, + hl/'dsl{N,,0) 
F3 = *6dJ(°)(iV„0), 
S = e*°°j(°)(A^c,0), 

F, = -e-*-(l + *io)d/i;^ A vol(4), (3.26) 
where dsg(A'c,0) is the deformed conifold metric and we defined 

POO 

= l^l^v / rfp' (sinh(4p') - 4p')"^^^ + hKs- (3.27) 
J p 

To understand the significance of this solution, notice that z/ 7^ leads to 

dF5 - /fs A F3 = 5 A S(4) ^ 0, (3.28) 

where the term on the right-hand side of this equation may be interpreted as the contribution 
from D3 brane sources smeared on the transverse directions. Indeed, the world- volume action 



for these sources arises in this limit, as discussed in subsection 2.2 
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Going back to the solution before taking the c — ?■ oo hmit, we can think of the term 



5 AS(4) in the second equation in (2.15) as a D3 brane charge density induced on the source 
D5 branes by the B field on their world-volume. Then we can compute the density of D3 
branes by integrating the B field pulled back to the world-volume cylinder wrapped by the 
D5 branes, where we put a cut-off at some radial distance. The result is also valid in the 
c — )■ cxD limit. Namely, we may define a running numbeinof source D3 branes as 



^ ' B (X NfC-^^ I e^'^+^'^dp', (3.29) 

cylinder 



where the factor of Nf comes from the overall factor in front of the action for the fiavour D5 
branes. Expanding this in the UV we get [qs = e*°°) 

Uf ~ ^,z/e^^/3 for p^oo. (3.30) 

The interpretation of this quantity becomes clear if we look at the asymptotic form of the 
warp factor in the standard radial coordinate r ~ e'^^^^. The leading term of the warp factor 
in the UV goes like 

z/r^ + A^^ loe r 

K - ^4 for r-^oo. (3.31) 

Expressing this in terms of the running number of source D3 branes nj, and running number 
of bulk D3 branes 



n ~ kNc ~ QsNllogr for r — )■ cxo, (3.32) 

this takes the form 

Qshy ~ -^-^ — for r — )■ oo. (3.33) 

This shows that there are precisely nf + n D3 branes in the background and reduces to the 
Klebanov-Strassler expression for u = 0. Notice that the running of the source and bulk D3 
branes is quite different, and in particular the former dominates the UV. We will discuss 
later the implications of this behaviour. In fact, this is the leading asymptotic behaviour of 



the warp factor (3.24) in the general case. 

The limiting solution is again singular near p = 0. However, this singularity comes entirely 
from the warp factor (3.27), while the metric dsg(A^c,0) is the smooth deformed conifold 
metric. This singularity is due to the fact that the D3 sources are distributed uniformly 
along the radial direction down to p = 0. 

^In the following formulas we will ignore numerical factors. 
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4 Comments on the field theory 



In this section we will discuss possible field theory interpretations of the solutions we pre- 



sented. As we will see, the solution of section 3.3 will be related to the mesonic branch of 



the Klebanov-Strassler theory, whereas a detailed field theory picture for the solutions with 
non-zero Nf and is more difficult to obtain. We find instructive to first recall the rela- 
tion between the gravity solutions and the field theory picture in the unfiavoured theories, 
as discussed in |5j. In that case the fivebrane solution contains one non-trivial parameter 
U oc c~^. Taking this to infinity corresponds to the near-brane limit, where a decoupled 
four dimensional field theory description exists. For any finite, non-zero value of U we can 
"rotate" the solution and add D3 branes and a new parameter p. In the new solution, the 
decoupling limit is obtained by taking /3 — t- oo. Thus U survives as a field theory parameter, 
and in particular it is proportional to the VEV of the U operator, partner of the conserved 
baryonic current. 

As already noted, the flavoured fivebrane solution has the same UV asymptotics as the 



unfiavoured one^ and is still charaterised by a parameter, denoted by c. The decoupling 
(near-brane) limit is obtained by sending c — )■ 0, where the asymptotic dilaton is linear, 
instead of constant. After the transformation that adds the D3 branes (for generic values of 
c), the warp factor will asymptote to a constant, unless again we take /3 — ?■ oo. We would 
like then to regard this as a limit in which a four dimensional field theory interpretation 
should be possible. To be more precise, let us write the leading expansion of the warp factor 
in the UV, after having dropped the constant term. Using the expansions in Appendix [B| 
we get 



h 



3cA^ 4p/3 , 1 

2 32 



{2N,-Nff{8p-1) 



297N] 



'Sp/3 



+ 0{e 



(4.34) 



and changing to the radial variable r ~ e^^^^ we have 



^ ^ 3e-^J- ^ ^r^ + 3(iVc-iV//2)Mogr | ^ 



(4.35) 



Thus in terms of the effective number of running D3 branes 



and 



n 



gs{N,-Nf/2f\ogr 



for r — > oo, 



(4.36) 
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Only the sub-leading behaviour is affected. 
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we see that agaiij"] 

gsh ~ '"'^ ' for r — )■ cxD. (4.37) 



rif + n 



Comparing with the hmiting solution discussed in section |3.3[ we see that the term propor- 
tional to n/ is due to the smeared source D3 branes, while the subleading term is due to 
the running bulk D3 branes. In particular, the former leading term does not give the usual 
asymptotic AdS with logarithmically running radius. This might suggest that a 4d field 
theory interpretation may not be available for these solutions. However, we can define the 
theory at some energy scale, i.e. at some radial distance r*. Once the field theory is defined, 
we may then follow its evolution in the far UV, as for the ordinary Klebanov-Strassler theory. 
Notice the present discussion is valid both for vanishing Nf or otherwise. 

The difference in asymptotics with respect to the usual logarithmic deviation from AdSs x 
T^'^ has an important implication in the dual field theory. In particular, we have that 
the would-be Goldstone boson fluctuation [29] around the solution of section 3.3 is not 



normalisable. This is easy to check by noticing that the ansatz in [29] is still valid in our 
case, and that the first integral in (3.27) of [29j is divergent with our asymptotic behaviour 
of the warp factor. Consequently, in our case we do not find a Goldstone boson associated 
with spontaneous breaking of the global baryonic symmetry U{1)b- Therefore U{1)b is not 
broken^, suggesting that our solutions always describe non-baryonic branches. The mode 
we are discussing corresponds to (infinitesimally) turning on the resolution parameter in 
the geometry [i.e. in the unwarped metric), which in our solution corresponds to turning 
on c~^. The fact that U{1)b is not broken strongly suggests that the operator U, defined 
as the partner of the conserved U{1)b current, does not have a VEV, hence we cannot 
associate the parameter in the gravity to this VEV. In fact, in the present situation this 
parameter may be interpreted simply as due to the back-reaction of the flavour D5 branes. 
In formulas, we have that = Nf/u, thus we see that introducing Nf D5 branes wrapped 
on a cylinder transverse to the S"^ of the resolved conifold turns on a finite size (as measured 
from infinity - see p]) for this two-sphere. More precisely, the parameter is proportional 
to the ratio of D5 to D3 sources. Notice that the D5 brane sources break the Z2 symmetry 
of the Klebanov-Strassler solution. We will elaborate on this point in section [5j 

Based on the above comments, we now turn to discuss more precise field theory duals to 



^^The factors of the dilaton combine with the factor in front of e^*dsg(-/Vt,, iVj) in (3.151. The factor of 
is canceled by an overall factor of c in the asympotic expression of ds§(-/Vc, Nf), so that the final metric 
does not depend on c at infinity. 

^^The absence of a Goldstone boson might be due to the fact that U{1)b is gauged. However this seems 
unlikely for our solutions. We thank Ofer Aharony for pointing this out to us. 
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the backgrounds we presented. We will start from the solution of section 3^ and then we 
will move on to the general solution. In this case, we will examine two distinct possibilities, 
pointing out pros and cons of both scenarios. 



4.1 Higgsing in the Klebanov-Strassler theory 



The solution ^'^ discussed in section 3.3 consists of an ordinary deformed conifold, with a warp 
factor comprising the standard source-less term Kks-, whose origin is the running number 
n of bulk D3 branes, plus a contribution arising from nj smeared source D3 branes. We 
therefore propose that the field theory dual to this solution is the Klebanov-Strassler theory 
with gauge group 
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SU{N^ + n + nf)x SU{n + nj) (4.38) 

on the mesonic branch. To understand this more precisely, let us first recall the structure 
of the moduli space of this theory, as discussed in detail in [27]. In the notation of this 
reference the generic Klebanov-Strassler quiver has gauge group 

SU{M{k + l) +p) X SU{kM + p). (4.39) 

The authors of [27] showed that the (quantum) moduli space of this theory consists of 
different branches, and may be written as a direct sum 

©^1 [Sym^a.o © Sym^-^-^a,! © Sym^-^^'^.s © ■ ■ ■ © Sym^^ Cr,k] (4.40) 

where p = kM+p. Each of the Cr^i denotes a deformed conifold, with deformation parameter 
depending on two integers r, /. The sum over r = 1, . . . , M is related to the M vacua arising 
in each of the terms, while the integer / labels branches of the moduli space arising at 
different steps of the cascade of Seiberg dualities. The power of the symmetric product is 
interpreted as the number of probe D3 branes moving on the corresponding deformed conifold 
geometry. The last term in the sum is the branch of lowest dimension, that is real dimension 
6p. However, if p = 0, the branch of lowest dimension is replaced by the baryonic branch, 
which is a copy of C [27] . 

Let us now consider a supergravity solution which incorporates the back-reaction of rij 
point-like D3 branes at distinct points on a warped deformed conifold geometry. This is 
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A closely related supergravity solution was discussed briefly in [T3] , although notice that the warp factor 



in our equation (3.27) differs from the corresponding one in equation (101) of 14J. Indeed the latter is a 
zero-mode of the source-less Laplace equation on the deformed conifold. 

^■^From now on we drop the tilde from Nc, Nf and denote them simply as N^, Nj. 
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different from the solution of section 13.31 The back-reaction of the source D3 branes is 
included by replacing the warp factor Hks hxs + J^jLi hj, where each hj is a solution to 
the Laplace equation on the deformed conifold, with a delta-function source. If the branes 
are placed at some point on the 5''^ at p = 0, we can interpret this solution in terms of 
a cascade of Seiberg dualities of the theory (4.39) down to the last step, where it becomes 
SU{M + p) X SU{p) and then it goes over to the smallest branch of (4.40), with p = Uf. 
However, if the D3 brane sources are placed at some finite radial distance, the natural 
holographic interpretation of the background is that the theory first undergoes a cascade of 
Seiberg dualities, and then it is Higgsed at an energy scale given by this distance. Then we 



are in some intermediate branch in (4.40). 



Let us now return to the solution of section 3.3, corresponding to smearing uniformly the 
Uf source D3 branes on the transverse geometry. Every time a D3 brane is crossed, the ranks 
of both gauge groups decrease by one unit. We may then interpret the gravity solution as a 
continuous process of Higgsing, occurring at all energy scales from the UV to the IR. Part of 
the change in the flux, in the geometry, and in the ranks in the field theory, is then due to 
this Higgsing, and is reflected in the u term in the warp factor. However, we still have the 
logarithmic part of the running, exactly as in Klebanov-Strassler, implying that the theory 
is also undergoing a cascade of Seiberg dualities. After resolving the singularity, it may be 
possible to stop the Higgsing at some finite distance, and proceed with the cascade. However, 
in our singular solution the Higgsing behaviour is dominating down to the IR, and therefore 
it seems that the Higgsed theory in the IR should be SU{Nc), where all of the D3 branes 
have disappeared. This possibility was entertained in [TS], who concluded that ultimately 
for the solution of |T3] the cascade interpretation is the correct one. On the other hand, in 
the presence of smeared sources like in our solutions, the Higgsing interpretation appears 
inevitable [15j . 

Standard supergravity computations, as for example the beta functions of the two gauge 
couplings, defined exactly as in [H], can be matched to the field theory computations, and 
are compatible with our interpretation. In summary, we have proposed that the gravity 



solution of section 3.3 corresponds to the mesonic branch of the Klebanov-Strassler theory 



with gauge groups given in (4.38), simultaneously undergoing a cascade of Seiberg dualities 
and Higgsing. 



^^A solution corresponding to coincident D3 branes localised at a particular point on the deformed conifold 
was discussed in ['28L 
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4.2 Higgsing in the Klebanov-Strassler theory with flavours 

Let us now address switching on a non-zero value for Nf. One way to think about the final 



background is as arising from the back-reaction of Nf D5 branee^^ placed in the resolved 
deformed conifold geometry of [Ij . However, K-symmetry requires to turn on a world- volume 
fiux on the D5 branes [SBj . Thus we have necessarily source D3 branes induced on the fiavour 
D5 branes. According to this point of view, the most direct interpretation of the solution 
with Nf and c^^ non-zero seems to be in terms of fiavour D5 branes on the baryonic branch 
of the SU{{k + l)Nc) x SU{kNc) Klebanov-Strassler theory. However, we cannot add Nf 
fiavour D5 branes without also adding extra (infinite) Uf D3 branes. On the other hand, 
from the point of view of the picture presented in section 4.1, turning on Nf and looks 



like a deformation of the mesonic branch of the SU{Nc + n + Uf) x SU {n + Uf) Klebanov- 



Strassler theory in (4.38). Below we will analyse two different scenarios for the dual field 



theory of our solution, based on these two different viewpoints, respectively. 
4.2.1 Explicit fiavour symmetry 

Using the fact that the Klebanov-Strassler theory is on the baryonic branch, it is natural 
to try to interpret the addition of infinitely extended D5 branes as true fiavours. In par- 
ticular, since the gauge coupling of the effective four dimensional theory on these branes is 



vanishing they should give rise to a global symmetry S'f/(A^/)fiavour- We will then consider 
a quiver of Klebanov-Strassler type, with an explicit global flavour symmetry group and 
quarks transforming in fundamental representations of the gauge groups, as well as of the 
flavour symmetry group. In particular, the presence of a single stack of flavour D5 branes, 
together with the global SU{2) x SU{2) symmetry of the background, suggest to consider 
the following gauge and flavour groups^ 



SU{N, + n + n/) X SU{n + Uf) x 5f/(iV//2)fl„, (4.41) 

and two pairs of quarks qi,qi, and q2,(j2, transforming under the two SU{2) factors. The 
total number of quarks is then Nf. Most of the discussion in this section applies with minor 
modiflcations to a general class of quivers, where the number of quarks of each gauge group 
may be different. In the following we will not specify a superpotential for the quarks. We 
will assume that there is a superpotential which allows to perform a sequence of Seiberg 
dualities. The theory we are discussing is depicted in Figure |6| where we denoted the ranks 



^^Reference [30] studied placing flavour D7 branes in the resolved deformed conifold geometry. 
^''This symmetry is typically broken to U{l)^f due to the smearing. 

^^We may consider replacing ^ N^. — Nf/2, without altering the main points of the discussion. 
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of the two gauge groups generically with A^^i and A^2- Note that this quiver is in the same 
class as the one discussed in [31] (cf. Figure 1 of this reference), where it was proposed to 
arise from a configuration with D7 branes. We study some aspects of this quiver in Appendix 

El 




Figure 6: A Klebanov-Strassler quiver, flavoured by the addition of Nj quarks. 



The first problematic issue that we face is that the ranks of the gauge groups in (4.41) 
include shifts by the Uf terms, arising from the source D3 branes induced on the flavour D5 
branes, as we discussed earlier. Then the resulting theory does not have a baryonic branch 
[27] , at least if the quarks do not have VEVs. Turning on VEVs for the quarks might change 
this. However, in our solution Nf is fixed, implying that the quarks cannot be Higgsed. 
In any case, we have not found a solution to the D-terms and F-terms of this theory, with 
or without VEVs for the quarks, that may be interpreted as a baryonic branch. Moreover, 
the absence of the Goldstone boson mode in the gravity, gives support to the fact that the 
theory is not on a baryonic branch. 

Building on the intuition gained from the discussion of the limiting case Nf — > 0, we 
would like to interpret the present solution in terms of a cascade of Seiberg dualities and 
Higgsing. Then the second problem that we face is that we are not able to match gravity 



computations to field theory computations in a natural way. The discussion around (4.35) 
suggests that the running number of bulk D3 branes is now given by n oc k{Nc — Nf/2), 
with the cascade step changing logarithmically in the UV as ~ {Nc — Nf/2) logr. If we 
assume that we can perform Seiberg dualities on the quiver in an alternate fashion on the 
nodes, it turns out that this is not the cascade pattern that the field theory follows. See 
Appendix [E] for details. Indeed, the difference in ranks, at each cascade step k in the quiver 
is Nc — kNf/2, thus it changes with the number of steps. This does not seem to be possible 
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in a string theory solution where the number of fivebranes is fixed, and only the number of 
D3 branes varies. Moreover, the rate of change of the ranks goes like k'^Nf, which is also 
not visible in our supergravity solution. One might speculate that there exists a particular 
pattern of Higgsing such that its contribution will exactly balance the terms of the cascade. 



in such a way that the final result detected by the gravity solution is that in (4.35). To us 
this does not appear to be natural, nor likely, and it is not clear how one could check such 
a conjecture. 

Finally, let us also mention that the supergravity beta functions obtained from identifying 
the dilaton with the sum of the gauge couplings, and the period of B with the difference, 
do not match with the field theory beta functions - see Appendix [E} In particular, in the 
field theory we have that the beta function for the difference of gauge couplings is 3Nc, 
independently of n,nf,Nf, which does not match the supergravity definition. Of course, 
there is no proof that these two calculations should agree, so this mismatch is not a rigorous 
objection to the validity of the present type of field theory scenario. We will discuss a little 
more beta functions in the following section, where we analyse a different field theory picture. 

In summary, if we follow the standard lore and add quarks transforming in the fundamental 
representation of the global fiavour symmetry, the resulting quiver field theory does not 
seem to be consistent with the properties of our gravity solutions. At least, if we work in 
the hypothesis that we can interpret everything in terms of cascade plus Higgsing. This 
interpretation is satisfactory in the limit Nj — > 0, when the quarks in the quiver of Figure |6] 
are absent. 



4.2.2 Emerging fiavour symmetry 

Assuming that a cascade of Seiberg dualties is taking place, we have seen that in the previous 
field theory picture, these are not refiected by the supergravity behaviour. In addition, the 
field theory beta functions do not match to the usual supergravity definitions, as one would 
have naively expected. It is then natural to wonder whether there exists a different picture, 
in which these features will match to the canonical expectations from the gravity side. We 
have been able to realise this only assuming that the field theory is a two-node quiver of the 
Klebanov-Strassler type, without exphcit fiavours added. In particular, let us then consider 
the following gauge groups 

SU{N, + n + nf)x SU{Nf/2 + n + nj). (4.42) 

After discussing how this proposal reproduces successfully various properties of the gravity 
solution, we will address what appears to be the main shortcoming of this picture, namely 
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the role of the global flavour symmetry. Notice that this theory reduces correctly to the 



field theory discussed in section 4.1 in the limit Nj — )■ 0, at fixed n/. Moreover, we can 
imagine removing all the n + nf D3 branes by undoing the "rotation" procedure. The theory 
then reduces essentially to that proposed in [71 [10] , modulo the fact that the second group 
appears to be gauged. We will return to this point shortly. 

First of all, the cascade of Seiberg dualities here matches to the gravity calculations 
in the usual way. We can define a running number n of bulk D3 branes by integrating 
A -F3 over the internal six-dimensional geometry (cut off at some distance in the UV), 
or equivalently from the corresponding UV term of the warp factor. Either way, one gets 
n oc gs{Nc — ^)^logr, and defining the cascade step k from the period of B over the 
two-sphere at large distance, we get 

k(xgs{N,- ^)\ogr n(xk{N,-^) (4.43) 



Then more precisely our Klebanov-Strassler quiver in (4.42) reads 



SU{N^ + k{N^ - Nf/2) + Uf) X SU{Nf/2 + k{N^ - Nf/2) + nj). {AAA) 

Of course the nj D3 branes may be still interpreted as being Higgsed, as we discussed. We 
now assume the validity of the following definitions for the gauge couplings and theta angles 

— + — = 7^e-^ — - — = 27re-'^ 1 - 60 , 
91 92 91 9i 

Bi + 62 = -27rCo, 61 - 62 = - / C2, (4.45) 



TT 



S2 



where we defined 



1 



B2 (4.46) 

S2 



and the integrals are performed on the two-cycle defined hy = [0 = 9 ^ip = 2Ti — (p^ip = -k]. 
It is well known that these definitions can be justified rigorously only in very special cases, 
but nevertheless they do capture the correct field theory quantities (at least in the UV) in 
a variety of cases. Evaluating the beta functions using the standard energy-radius relation 
we then obtain (in the UV) 

A = -/32 = 3(iV, - ^), AGi = -AOs = a{2N, - Nf) . (4.47) 
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Assigning (in the UV) anomalous dimensions and R-charges to the bi-fundamentals, exactly 
as done in [H] 

lA^ =lB, = -\, Ra= Rb, = \ (4.48) 

we find that the field theory calculations reproduce the result in ( |4.47 ). We note that the 
computation of the holographic central charge given in Appendix [F] is perfectly compatible 
with the quiver theory that we are discussing here. 

The remaining issue to address is the fate of the expected SU{Nf/2) global fiavour sym- 
metry in this theory. We would like to propose that we can imagine the quarks arising in 
the IR, after the cascade of Seiberg dualities and the Higgsing have taken place. In other 
words, these should be the bi-fundamentals Ai,Bi, in a configuration in which one gauge 
group is very weakly coupled, and can be thought of effectively as a flavour. More precisely, 
writing the matrix indeces explicitly (suppressing the SU{2) x SU{2) indeces), we can deflne 
a^j = APj, Vp = B\ and q^a = A^ai Q^p = B°'p, where p labels Nc + n -|- nf , j labels n + nf 
and a labels Nf/2. At every cascade step Nc and Nf stay flxed while n + Uf decreases, and 
similarly, at every Higgsing step the bi-fundamentals a, b get mesonic type VEVs. Hence this 
splitting is well deflned at any step. At the end of this process, the a, b bi-fundamentals will 
disappear, leaving the q, q "quarks" . This is similar to the discussion in Klebanov-Strassler 
(see around eq.(109) in [S]), where the authors proposed that the addition of a small num- 
ber p of probe D3 branes at the bottom of the cascade should correspond to a SU{M + p) 
gauge theory, with 2p "flavours". Indeed, if we take A^//2 ^ Nc, at the bottom of the 
cascade/Higgsing process, the gauge coupling of the SU{Nf/2) group will be much smaller 
than that of the SU{Nc). Thus, at least in this case, this proposal is reasonable. 

Let us make a few more comments in support of this idea. First of all recall that the global 
flavour symmetry in the dual theory arises from the gauge symmetry on the world-volume of 
the flavour D5 branes. In particular, global currents are given by fluctuations of the gauge 
flelds on the branes. However, if there were no normalisable fluctuations, then we would not 
have any state charged under this symmetry, and the global symmetry would be effectively 



gaugec ^ Given our asymptotics, it may be the case that all fluctuations of the gauge flelds 
on the flavour branes are non-normalisable. It would be interesting to prove or disprove this 
statement. We also note that in our solution the dilaton e~* is divergent at the origin (see 
(3.22)), suggesting one of the gauge couplings vanishes in the IR. We do not know whether 

^^Notice that before the "rotation" the Nc D5-Nf D5 sohition is not decoupled from gravity, hence any 
putative global symmetry should be gauged. Only in the c ^ limit does the sohition correspond to a field 
theory with Nf flavours. 
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we can trust this, or if it is an artifact of the IR singularity. Finally, recall that in the 
ordinary cascade one can get close to points where one of the couplings vanishes - see [3^ - 
and we can think that the Higgsing may eventually freeze the running at this point. 



Interestingly, note that in the case Nf = 2N(. the quiver (4.44) becomes actually a 
Klebanov-Witten [33] quiver 

SU{N^ + Uf) X SU{N^ + Uf). (4.49) 

It is satisfying that in this case the gravity calculations show that the beta functions and 
the anomalies vanish in the UV. Moreover, the cascading behaviour disappears from the 
solution, as can be seen for example from the leading behaviour of the warp factor. This 
solution may be then describing a Klebanov-Witten theory, where conformal invariance is 
broken by the Higgsing effect, represented by the leading term in the warp factor. As can be 
seen neatly from the presentation in Appendix |Dj in this case a Z2 symmetry is restored in 
the solution. The U{1)r (corresponding to ■?/'—)■?/' + 5) is broken to Z2 but, differently from 
[TT| [13], it is not broken in the UV, as expected from the fact that this is not anomalous in 
this case. 

We note that there exists also an analytic solution, which does not correspond to setting 
Nf = 2Nc in the solution discussed in the rest of this paper. In this solution the ?7(1)_r is 
preserved, while there is no Z2 symmetry. We briefly discuss this solution in Appendix [C] 
Solutions with Nf = 2Nc and their possible desingularisations are clearly interesting in their 
own right. We leave further studies of these solutions for future work. 



The two scenarios illustrated above are by no means exhaustive and perhaps a more 
sophisticated set-up is necessary to obtain a better understanding of the field theory dual to 
the background we have presented (if indeed this exists). 



5 The Z2 in gravity and in field theory 

In this section we will discuss a Z2 symmetry in gravity, as well as its field theory interpre- 
tation. 



5.1 Back to gravity 

Following |5] , we can define an effective (running) resolution parameter a^ff as the size of the 
two-sphere measured from infinity. In the unwarped metric this gives ~ [Nc — Nf/2)p. 
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See Appendix D In the warped metric after the rotation, this becomes multiphed by the 
parameter c~^, hence we can write this as 



Nf 



(5.50) 



We can then imagine starting from the solution of section 3.3 and increase Nj from zero 



keeping u fixed. Interestingly, the effective resolution will then increase, reaching a maximum 
at Nf = Nc and then will vanish again at Nj = 2N(.. Increasing Nf further, the resolution 
parameter will become negative. There is another solution, obtained by a Z2 symmetry, 
which in Figure [7] corresponds to reflection along the horizontal axis. The points Nf = 
and Nf = 2Nc are special, as we now discuss. 





-60 L 



Fi gure 7: Plot of the effective resolution parameter Oi^ff as a function of Nf at fixed u and 
(large) p. The two branches are exchanged by a Z2 reflection. At Nf = the Z2 symmetry 
X is restored. At Nf = 2Nc the reflection symmetry TZ is restored. 

In the gravity side there is the Z2 action which interchanges the two-spheres defined as 



TZ : 9 ^ 9, ip ^ ip. 



(5.51) 



This can be combined with the change of sign of the three-forms — )■ — i/3, -F3 — > — F3 (the 
center of SL{2,Z)) to define the following Z2 action 



(5.52) 



Notice that the effective resolution parameter a^ff changes sign under TZ and X. These 
operations have been discussed in the context of the Klebanov-Witten [33] and Klebanov- 
Strassler [22] theories. Let us briefly summarize the behaviour of various gravity solutions 
under the X (and 71) action: 
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• Unflavoured. Unrotated. For any value of c < oo the effective resolution is non-zero, 
hence the metric is not invariant, while the F3 is invariant. Therefore the full solution 
is not invariant under X (neither invariant under TZ). 

• Unflavoured. Rotated. For any value of c < cxd the effective resolution is non-zero, 
hence the metric is not invariant, while the F3 is invariant. Therefore the full solution 
is not invariant under X (and TZ). For c = 00 we have the Klebanov-Strassler solu- 
tion: the effective resolution vanishes hence the metric becomes invariant. The F3 is 
invariant. Hence the full solution is invariant under X (but not under TZ). 

• Flavoured. Unrotated. For any value of c and Nf 7^ 2Nc the effective resolution is 
non-zero, hence the metric is not invariant. The F3 is not invariant. Hence the full 
solution is not invariant under X (and TZ). If we exclude the solutions in Appendix [C| 
for any value of c and Nf = 2Nc the effective resolution vanishes hence the metric is 
invariant. The F3 is not invariant under X but it is invariant under TZ. Hence the full 
solution is not invariant under X but it is invariant under TZ. 

• Flavoured. Rotated. For any value of c < 00 and Nf 7^ 2Nc the effective resolution is 
non-zero, hence the metric is not, and the F3 is not invariant. Hence the full solution 
is not invariant under X (and TZ). For any value of c and Nf = 2Nc the metric is 
invariant. The F3 is not invariant is under X but it is invariant under TZ. Hence the 
full solution is not invariant under X but it is invariant under TZ. For c — )■ 00 and 
Nf — the effective resolution vanishes hence the metric is invariant, and the F3 is 
invariant. Hence the full solution is invariant under X (but not under TZ). 

5.2 Field theory 

In the Klebanov-Strassler field theory it is argued in the literature that X acts as follows: 

X : Ai^ B}, Bi^ A\ (5.53) 

which is interpreted as simultaneous exchange of Ai and accompanied by charge conju- 
gation (without exchange of gauge groups, which does not make sense for unequal ranks). 
Recalling that the Klebanov-Strassler superpotential y&W = Tr\AxBxA2B'i\ — Tr\AxB2A2B^ 
we see that 

X: W^W*. (5.54) 
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Since we have to take the real part of W in the Lagrangian, we see that this is invariant 
under X. On the other hand, we have 

X : K = Y^ Tr[AiA\ - BjBi] -U, (5.55) 

i 

which explains why the Klebanov-Strassler point does not break X, while for generic points 
on the baryonic branch, this is broken. 

We now want to study the action of X when there are quarks, coming from flavour branes. 
Notice that under X the smearing form S(4) changes sign. This means that the flavour 
fivebranes change orientation, and suggests that we should still apply charge conjugation on 
the quarks. The interchange of the bi-fundamentals Ai and Bi under X is dictated by the 
fact that the interchange of the two-spheres is a "geometrical" operation [33J. To see this. 



it is convenient to write the (resolved) conifold equatioE^° as ^^+^2+^3+^4 = 0. By 
writing the Wi coordinates explicitly in terms of angular coordinates it is then easy to verify 
that [33] 

X : Wi ^ Wi, W2 — )■ W2, Ws — )■ W3, W4 — —W4. (5.56) 

Then a further change to the variables Zi defined as Wi = zi + Z2,iw2 = Zi — Z2,W3 = 
zs + Z4^, iwi = z^ — Z4, where zi = AiBi,Z2 = A2B2, Z3 = A1B2, z^ = A2B1, makes clear why 
X interchanges Ai with Bi. On the other hand, the quark fields come from the open string 
stretching between the flavour fivebranes, and they are not related to the geometry. It is 
therefore plausible that the quarks should not be interchanged under X. Then generically, if 
we imagine that our field theory has bi-fundamentals a^, bi and quarks qi,qi, we can postulate 
the following action 

X : tti^ b], bi a|, Qi qj, Qi qj. (5.57) 

In order to see whether this is a symmetry of a field theory, we need to specify the 
Lagrangian and the vacuum. Let us concentrate in the field theory scenario discussed in 



4.2.2 As explained earlier, we are thinking that the "quarks" are part of the bi-fundamental 



fields. In matrix notation, the splitting discussed in section 4.2.2 can be written as 

A = {a g), B^ = [b^ f), (5.58) 

where a, b^ are (Nc + n + Uf) x (n + Uf) matrices and q, q^ are (Nc + n + nf) x Nf/2 matrices. 
When we write these expressions we are making manifest only a SU{Nf/2) x SU{n + nj) 

^°This calculation is justified by thinking about the brane set up before the back-reaction. 
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sub-group of the SU{Nf/2 + n + nj) gauge symmetry. If we apply SU{Nf/2 + n + nj) 
rotations, the a and q, and b and q in ( |5.58 ) will mix. However, it is legitimate to define 



how a symmetry (namely X) acts in a convenient gauge. Notice also that in our picture by 
definition the a, b will get mesonic-type VEVs during the Higgsing process, while the q, q 
will not get VEVs. This gauge choice is then also convenient for this purpose, since as we 
cascade/Higgs, the rank of SU{n + nf) will decrease, while SU{Nf/2) will stay fixed. 
Writing the Klebanov-Strassler quartic superpotential in terms of a, b, q, q we have: 



W = Wi + W2 + W3 (5.59) 



where 



Wi = Tr[aibia2b2] - Tr[ai62a2&i] 
W2 = Tr[gi 51^2^2] - Tr[qiq2q2qi] 

W3 = Tr[qiqia2b2] + Tr[aibiq2q2] - Tr[qiq2a2bi] - Tr[aib2q2qi]- (5.60) 
Then we have (the real part is understood) 

X : Wi^ Wi, W2 -W2 (5.61) 

and 

X : Pl/g = Tr[gigia2&2] + Tr[aibiq2q2] - Tr[gig2a2&i] - Tr[aib2q2qi]- (5.62) 

Notice that the transformed is not related simply to the initial W3. However, this still 
makes sense, since in the transformed theory we reversed the arrows on the quarks. 

We can ask whether the mesonic VEVs turned on in the Higgsing process can break the 
Z2. In fact, if we think of adding D3 branes at points in the deformed conifold, in general the 
resulting solution will break all the symmetries of the deformed conifold. However, since in 
our set up the source D3 branes are smeared uniformly, all the symmetries that would have 
been broken at a generic point on the mesonic branch, are restored. This is obvious from the 
fact that our limiting Klebanov-Strassler-like solution has all the symmetry of the deformed 
conifold, including the Z2. Hence in our solution we should think of the Z2 breaking as due 
only to the quarks, and not to the mesonic VEVs. 

As we discussed, in the case Nf = 2Nc the metric and the three-forms become invariant 
under the interchange of two-spheres, hence a different Z2 symmetry appears. Since now 
the ranks of the gauge groups are equal, it seems natural to assume that this Z2 action in 
the field theory corresponds to the interchange of Ai and Bi, together with the exchange of 
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the two gauge groups, as discussed in [33JqJ In our case this operation exchanges the role 
of "flavours" and "colours", which of course can make sense only in the quiver (4.49). 



6 Discussion 

In this paper we have used a solution generating transformation, applicable to a large class 
of supersymmetric Type IIB backgrounds, to construct a new family of solutions generalising 
the resolved deformed conifold of Butti et al [1]. Using this method, we can take any solution 
to the torsional superstring equations [31 [191 Elj and generate a solution where various RR 
and NS fields are turned on. The method may be applied to solutions which include the back- 
reaction of smeared source branes, usually referred to as fiavour branes. In particular, we 
have applied this procedure to a solution representing a system of Nc D5 branes wrapped on 
the two-sphere inside the resolved conifold, with addition of Nf flavour D5 branes, wrapped 
on a transverse inflnitely extended cylinder. The flnal solution is then a warped resolved 
deformed conifold, modifled by the back-reaction of the extra flavour branes. 

The flavoured solution differs qualitatively from the unflavoured one in two ways. Firstly, 
it is singular in the IR. Secondly, the UV asymptotics is not the (logarithmic) Klebanov- 
Strassler one. Although we have not addressed the resolution of the IR singularity in this 
paper, we expect it may be resolved by considering a proflle for the flavour D5 branes that 
vanishes smoothly in the IR. The different behaviour in the UV is induced by the presence of a 
uniform distribution of source D3 branes, smeared on the transverse geometry, up to inflnity. 
We have explained that these D3 branes are induced by the presence of the flavour D5 branes 
in a geometry with a non-trivial i?-fleld. In other words, the "rotation" procedure has the 
effect of adding bulk D3 branes, coming from the original colour D5 branes, and smeared 
source D3 branes, coming from the original flavour D5 branes. 

This geometric set-up leads to an unusual picture in the dual fleld theory. First of all, the 
different asymptotic behaviour of the warp factor implies that the zero mode associated with 
changing the parameter is here non-normalisable. This suggests that our gravity solu- 
tions correspond to vacua where the global baryonic U{1)b is not spontaneously broken. In 
particular, assuming that the fleld theories are some modiflcations of the Klebanov-Strassler 
theory, we have proposed that our solutions correspond to mesonic branches. This matches 
nicely with the fact that we always have n f extra D3 branes in the geometry, which are run- 
ning because of a Higgsing effect. We can then explain the parameter purely as arising 

^-"^This operation changes the overall sign of the superpotential. This may be reabsorbed by the symmetry 
T, which by definition leaves invariant the lowest components of the superfields Ai, Bi 133j. 
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from the back-reaction of the Nf flavour D5 branes in the geometry. In particular, since 
changing c"^ corresponds to changining Nj (at fixed n/), this is a change of theory. This 
is rather different from the interpretation in the case without any source D3 or D5 branes, 
where this parameter is related to the classical VEV of the operator W, which is the partner 
of the baryonic U{1)b current. This different interpretation of a background in the presence 
of extra sources should not be too surprising. Indeed, the Klebanov-Strassler geometry, plus 
some probe D3 branes is indeed dual to the field theory in the mesonic branch. On the 
mesonic branch we do not expect a fuzzy two-sphere emerging from the field theory. The 
asymptotic growth of the warp factor in our solution indicates that the effective radius of the 
two-sphere in the closed string metric is never much smaller than the B field suggesting 
that a fuzzy sphere interpretation is in this case not applicable [M]. However, we have not 
checked this in detail. 

The general picture that we have proposed is that in our solutions the decrease of five-form 
flux, and of ranks in the gauge theories, is due to a sequence of steps involving both Seiberg 
dualities and Higgsing. In the case where we turn off the Nf flavour D5 branes (but not 
the n/ D3 sources) this picture is robust. Turning on Nf led us to discuss slightly more 
exotic scenarios. The standard lore about "flavouring" fleld theories instructs us to pick 
an unflavoured theory and add quarks transforming under a new global flavour symmetry. 
However, we have seen that in the present case this recipe leads to a fleld theory that 
fails to reproduce the main properties of the gravity side. In particular, it does not seem 
possible to reproduce the pattern of Seiberg dualities. This is perhaps not very surprising, 
given that flavoured versions of the Klebanov-Strassler theory of the type we considered in 
section 4.2.1 actually were shown to arise from constructions involving D7 branes [31]. We 
have then analysed the possibility that the flavour symmetry might actually emerge only 
in the IR. Although this may be an unorthodox suggestion, we have explained that such 
proposal matches all the expected gravity predictions. However, checking this statement in 
our solution is difficult, since this is singular in the IR. 

We have indicated that the singularity in the IR may be resolved. Hopefully further studies 
of non-singular versions of our solutions will clarify the validity of the ideas presented in this 
paper. 
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Appendices 

A Type IIB equations of motion with sources 



In this appendix we state the equations of motion coming from the action (2.13). Let us 
first rewrite the total action, comprising the Type IIB supergravity action plus the super- 
symmetric source action: 

S = SjjB + •S'sources (^-l) 

where 



SiiB = I V^[R- ld^<^d^<^^ C(4) A F(3) A H 



Iff 1 ^ ^^-'^ 

" 2 / ( '''^^(i) ^ + ""'"^ ^ *^ + A *F(3) + -F(5) A *F(5) 

and 

^sources = ' J (e^^+*/\ol(4) A ( COS (c^^ J + siu (e-^/^B) - C(6) + ^(4) A A S(4) (A.3) 

In the following we will set = 0. The modified Bianchi identities for the fluxes read 

dH = 

dF(3) = S(4) (A.4) 
dF(5) = A F(3) + S A S(4) 
The equations of motion for the fluxes read 

d (e-* * H) = F(3) A F(5) + sin Ce'^vol(4) A S(4) 
d (e* * F(3)) =-HA F(5) 

Notice that in the equation of motion for H the term coming from C(4) ABA S(4) in the 
source action is exactly cancelled by a contribution from the Chern-Simons term of the Type 
IIB supergravity equations. We then define the following notation 

^ipy\p) = ^^""■•"''A^,,..^^ (A.6) 

We also have that 

j uji^p) A A(io-p) = - j v^Aj(*a;) (A.7) 
Using these, we can write the dilaton equation of motion as 

d, [V^gg^^d,^) = -^e*F(|) - ^e-^H' - le^/^.^. * (cos Ce^^vol(4) A J) (A.8) 



-g - - ' 12" '(3) 12" 2 
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Finally, the Einstein equation is 

J ^ P P P1P2P3P4 

96 (A.9) 

- ^e^^+*/2 cosC (2S,,,,,,3 * (vol(4) A J)/^"^''^ - 3(7,.S(4). * (vol(4) A j)) 

- ^ sinCe^^ ((5 A H(4));.p,...p, (*vol(4));"-''^ - 60(7^,(5 A S(4))j (*vo1(4))) 

B More on the solutions 



In this appendix we discuss some properties of the master equation (3.7) and we present an 
analysis of the asymptotic form of the solutions discussed in the main body of the paper. 
Moreover, we desribe a systematic method for constructing these solutions recursively. Let 



us start by noticing that the master equation (3.7) can be written in the compact form 



4s (P' + Nf] 



where 



Defining 



we obtain 



(B.2) 



sinh2(2(p-p,))' 
P = 2 y dps-^P, (B.3) 

VsdpisdpV) - ^{sdpVf + 2NfsV + Q2 = 0. (B.4) 



Moreover, introducing the variable 

t^2 I dps-' = I ^ (B.5) 

J [|(sinh(4(p-po)) -4(p-po)), p^ > -oo, 

we have P = V and the master equation takes the form 

AVV -V^ + giit)V + goit) = 0, (B.6) 

where the dots denote derivatives with respect to t and the functions gi and g2 are given by 

g^{t) = 2NfS, g^it) = Q\ (B.7) 
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Equation (B.6) is equivalent to (3.7), but it is only quadratic in V, which considerably 
simplifies the systematic analysis of the equation. On the other hand, for po > — oo the 
functions gi{t) and go(t) are only parametrically known in terms of p, but this will not be a 



serious obstacle in the following analysis. Notice also that (B.6 ) is manifestly a generalisation 
of eq. (A. 6) of [5J to the case with Nf ^ 0, upon identifying V = N^f and t = t. The 



simplified master equation (B.6) can be derived from the simple action 



S = j dtV-^'^ {v^ -gi{t)V + 



9o{t) 



(B.8) 



which again is a generalisation of the Lagrangian (A. 7) of [5J for Nf ^ 0. The corresponding 
Hamiltonian is given by 



H = V-'/'{^Vn^j, + giit)V-go{t) 



where the canonical momentum conjugate to V is 

It follows that the solutions can be obtained from the Hamilton- Jacobi equation 



(B.9) 



(B.IO) 



with 



vr-p 



dS_ 
dV' 



(B.ll) 



(B.12) 



Obtaining a solution S{V,t) of (B.ll) leads to a first order equation for V{t) via the iden- 
tification 



2-p-l/2p 



95 
dV' 



(B.13) 



B.l Recursive construction of the solutions 



It is known [3110112] that the solutions of (3.7) fall into two possible categories according to 



their UV behaviour: they either lead (via (3.6)) to an asymptotically linear or asymptotically 



constant dilaton. Since the "rotation" formulas (3.15) require that the dilaton is bounded 



from above, it is only the solutions with the asymptotically constant dilaton that can be 
rotated. As we have seen, these solutions behave as P ~ ce^''/^ in the UV, for some strictly 
positive constant c. In this section we present a systematic way of constructing all solutions 
with these UV asymptotics in an expansion for large c [9]. 
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Instead of constructing these solutions by expanding P, as was done in [9], we will use 



(B.6) and write 

V{t) = dp{t), (B.14) 
where c = 3c/2, has been chosen such that (for po = 0) 

V{t) ~ ^ P(p) ~ ce^^'/l (B.15) 

Expanding 

oo 

p(t) = ^c->„(t), (B.16) 

n=0 



and using (B.6) we get the sequence of equations 



4poPo-p^ = 0, (B.17) 

^ / n-l \ 

2poPn - PoPn + 2poPn = ~2 ( ^1^"-^ + ^O^'^^ + ^1 (4PmPn-m - VmPn-rn) 1 , n > 0. 

The general solution of the equation for po is 



Po = a 



{t-t,fl\ (B.18) 



where a and to are constants, a can be absorbed in c and we choose to = so that t — )■ as 
p — )• Po- Using then po = ^^''^ in the rest of the equations, we find that the general solution 
is given by 

p^{t) = t^l^ f dt't'-'/^Rn{t') - f dt't'-'/'Rn{t'), n>0, (B.19) 
Jt± Jt~ 



where 



^ / n-l \ 

Rn{t) = -^t^^^ I gi{t)Pn^l + go{t)Sn2 + {^PmPn-m " PmPn-m) j , (B.20) 

\ m=l / 

and are arbitrary constants. 
Some comments are in order here: 

• Given this expansion for V{t), P{t) can be obtained by term-wise differentiation as 

oo 

P = p = ^gi-"p„. (B.21) 

71=0 
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By setting t+ = oo one ensures that p„ = o{t^^^) as t — )■ oo, i.e. that p„ for n > is 
asymptotically subleading in the UV compared to po. In fact, one can prove a stronger 



statement, namely that the expansion (B.16) is also an asymptotic expansion for the 



solutions that asymptote to the deformed conifold as t — t- oo. In particular, with the 
choice t+ = oo and using the fact that 

1 



g,{t) = ^{2N,- Nff log' {8t) + 



one can show iteratively that 

Pn{t) = O (t(4-")/3(iogt)'=W) , as t ^ OO, 



(B.22) 



fB.23l 



where k{n) is a positive integer. Hence, Pn+i is asymptotically subleading relative to 
Pn, i.e. not just compared to po. 



In the flavour-less case the expansion (B.16) is uniformly valid for all t. In particular 



with a suitable choice for t~ it correctly reproduces the IR asymptotic behaviour of 
the Butti et al. solution [1], as well as the UV behaviour. As we demonstrate below, 
this allows us to determine the constant that parameterises the IR asymptotics in an 
expansion in large c, the parameter that governs the UV asymtptotics. 

Although the integrals involved in the expressions for p„(t), n > 0, cannot be evaluated 



explicitly, the expansion (B.16) is converging rather fast as we will see in the plots 
below. 



Inserting the expansion (|B.16|) in the expression for the dilaton in (3.6) we obtain 

3 



-2<S> 



(l + l^tit-'/'p^] 



2c2 



5/3 „2 



where 



-2*oo ^ -2*,~3/2 

9 



(B.24) 



(B.25) 



This formula is useful for evaluating the warp factor h. 



B.2 The unflavoured solution 



Let us now examine in more detail the expansion (B.16) in the case Nf = 0. Firstly, since 
gi(t) = in this case, it follows by induction that all odd terms P2k+iit) vanish. Moreover, 
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the unique value of tgfc for which the expansion (B.16) leads to a regular solution as t — )■ 
is tgfc = 0- With this choice of the integration constants the expansion (B.16) for the 



flavour-less case is uniformly valid for all values of t G and so one can immediately extract 
both the UV and IR behaviour of the solution. 

UV asymptotics 



As t — )■ oo we have 

9o{t) 



log't 



^(log2(8t)-41og(8t)+4 + ^ 



Using the recursion relations (B.19) then one computes 



Differentiating this with respect to t and using the relation 



(B.26) 



P2it) = ^t'/' (^log^(8t) - 71og(8t) + y) + Oit-'/Hog'it)). (B.27) 



t 



p + 0{e ^^), as p oo, 



we obtain the following asymptotic expansion for P{p) 0: 

c \ 16 / 



Via (3.5) then we obtain 



,4(<I.-4.„) 

a 



ce'^p/s + jVe(l - 2p) + 0{e'^^/^), 
3 6c 



2e-2^ - 



2N, 



[1 - 2p)e-i°^/3 + 0{e 



(B.28) 



(B.29) 



(B.30) 



IR asymptotics 



As t — )■ we have 



9oit) = ^(i3tr/'-li3tr+oit'/' 



(B.31) 
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With the choice t2^. = discussed above then, it is easy to show that p2k{t) admits an 
expansion of the form 

oo 

P2.(t) = $^p2t^, (B.32) 

£=0 



where 



/■oo 

= 1, = 0, £ > 0, pf^ = - / dt't'-'/^R2kit'), k>0. (B.33) 

Jo 

Hence, 

oo oo oo 



fc=0 £=0 £=0 



In principle, one can obtain this expansion directly from the recursion relations (B.16) as we 
did for the UV asymptotics, but each term in the IR expansion will contain an infinite sum 
over powers of c: 

oo 

= (B.35) 

fc=0 

Crucially, all these sums can be expressed algebraically in terms of the leading term hi = p^^\ 
One then obtains the expansion 

q/V2Q-2/3 9iV2 f81Ar23-4/3 - I60hj) 

V = hit^l^ - e ^ . ^t«/3 ^ 0(^10/3)^ gg) 

40/ii 44800/if 



Alternatively, this expansion can be obtained by inserting (B.34) in (B.6). Differentiating 
this with respect to t now and using 

we obtain the IR expansion [3, [9j 



where we have defined 

34/3 

hi = —hi. (B.39) 
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Using this expansion then we get via (3.5) 



=2g 



=29 



+ 



4 
45 



,4($-<I.o) 



hi 4 

\ 

2 15 

+0{p% 
hi 2 {hi 
2 5hi 
32 



6hi + 15Nc - 
2N? 



~hr 



3hi - 5Nc 



p' + 0{p 
]{3hf 



rohfNc 



lUhjN^ - 32N^) p^ 



+ 



hi 

8 {3hi - AhjN'^ 



1575/1? 



1 



1+ -2 + 



64iV2p2 i28iV2 
3hi 



315/if 
-15hl + 124A/'2)p4 



+ Oip' 



+ 0{p' 



P' + 



405hf 

2 {75hl - 232hlN^ + imhiN^ + UN^ 



Dp' 



Abhl 



(B.40) 
+ 0{p'). 



Before we move on to the flavoured solution, let us point out a couple of interesting 
properties of the unflavoured solution, which will not be shared with its flavoured version. 



Firstly, the fact that the expansion (B.16) is uniformly valid for all t, implies that not 



only both the UV and IR asymptotic expansions can be directly obtained from the large-c 
expansion as we discussed above, but we can also relate the IR parameter hi with the UV 
parameter c. Namely, we have 



hi 



34/3 



E 

fc=i 



3c 
2" 



l-2k 



dtr'"'^R2k{t) 



(B.41) 



which provides a systematic way to obtain hi in an expansion for large c to any desired order. 
Since the dilaton goes to a constant both in the UV and in the IR, as seen in the above 



asymptotic expansions, we can view (B.41) as a relation between the UV and IR values of 
the dilaton. In fact, for Nf = 0, it follows from (3.6) and (B.l) that the Hamiltonian (B.9) 
is essentially the dilaton [5]: 



i7 = 4e-2(*-*°), 



(B.42) 



and hence (B.41) relates the IR and UV values of the Hamiltonian. As can be seen in Fig. 



8 the expansion (B.16) and the resulting relation (B.41) converge very fast. 



B.3 The flavoured solution 



Let us now try to construct the flavoured solution of (B.6) that reduces to the unflavoured 
solution presented above as Nf — 0. As we shall see, in the flavoured case the expansion 



(B.16) is not as useful as in the unflavoured case since it breaks down in the IR and hence it 
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Figure 8: The function and the dilaton plotted using the first three orders in the 



expansion (B.16) (0(c) red, (9(c~^) blue, and (9(c~^) green) are compared to the numerical 
solution. The plots are for the values c = 79.370, Nc = 10 and hi = 223.3, the last of which 
is determined numerically. Note that the numerical v alue fo r hi is very close to the leading 
order approximation hi ^ 40/3^^^ in the expansion (B.41) relation, which, as seen in the 
plot of the dilaton, converges very fast. 



is not uniformly valid for all t. It nevertheless provides a useful representation of the solution 
in the UV. 

UV asymptotics 



The derivation of the UV asymptotics from (B.16) proceeds as in the unflavoured case. 
As t — )■ oo we have 



goit) = ^(2iV, - Nff (^log2(8t) - 41og(8t) + 4 + O 



(B.43) 



Using the recursion relations (B.19), we compute 



9A^f 



(B.44) 



P2(t) = It'^' (^^(2iV, - Nff (log^(8t) - 71og(8t) + ^) 



9A^ 



O 



Differentiating these with respect to t we now obtain the asymptotic expansion [U] : 

8iivr 



64 



+0{pe 



(B.45) 
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Via (3.5) these lead to 



.2g 



,29 



,4($-$„) 



+ - (2iV,-iV^)(2p-l) + 



+ 0(e- 



-4p/3N 



ce4p/3 + 1 ^(2iv, - iV^)(l - 2p) + + 0{e-'^") 



3 2 24c 

2c3 V c ' 16c2 



(2iV,-iV;)^(4p-5)^ 



9iV 



-4p/3 



-8p/3N 



+ 



((2Ar, - Ar^)2(l - 8p) + 297A^J) e"^^/^ + C(e" 



'4p^ 



26-2"- -(2A/'e- A^/)(l -2p)f 



'lOp/3 



+ 0{e 



-14p/3N 



(B.46) 



/i? asymptotics 



Turning to the IR asymptotics, we need the expansions 



2Nr-N 



9o{t) 



g^{t) = 2Nj [ {3t) 



9 \ 5 



-2/3 



_ 1 + 0(t2/3)^ _ 



(B.47) 



Inserting these in the recursion relations (B.16) we find for small t 



Nf 



32/38/11 3V332/i2 



iV2 \ jV3 

^ ^ log t + ^t^/^ log^ t + . 



576/if 



(B.48) 



It can be shown that this expansion is what one obtains by first solving (B.6) perturbatively 
in Nf around the unfiavoured solution of the previous section and then looking at the IR 
behaviour of this solution. However, it is obvious that this expansion is not valid all the way 
down to t = as it breaks down at t ~ exp(— 32/'^2/ii/A^j). It follows that the expansion 
(B.16) breaks down in the IR for Nf ^ and (B.48) is not the correct IR behaviour. To 



recover the correct IR behaviour of the fiavoured solution we need to go back to (B.6) and 
look for an asymptotic expansion as t — > non-perturbatively in Nf. 



It can be shown that such an asymptotic solution of (B.6) is given by 

32/32 r 



V{t) 



Nf 
327^ 



^4/3 



logt + 



, 1 , , , 1 log(-logt) 

— hi - - log - log t - — — 

Nf 4 16 logt 



1 ^32/32- 5 
-hi - 



4 V N 



f 



1 



A logt 



O 



log(-logt) 

(l0gt)2 



(B.49) 
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where again hi is an arbitrary constant. Differentiating witli respect to t and using 

8p3 32p^ 256p^ 



t = — H '— + 

3 15 315 



(B.50) 



we obtain 



wliicli leads 


via 


(3.6) to 


e^^ = 




g logp + /ii 


e^^^ = 


ll 


logp + /li 






logp + /li 


g4(<I.-4.o) ^ 




27 




2Ar|(-logp)3 



AT. 



9 



log(-logp) + 



log(-logp) + 



log(-logp) 

logp 
log(-logp) 



_Ni_Ni 
3 9 
log(-logp) 



log(-logp) + 



logp 

log(-logp) 



logp 



+ 0(p3logp), 

(B.51) 

p2 + 0(p3logp), 
+ C(plogp), 

+ C>(p2logp), 



i-2pMi. (2^-^^ 



4 log p \ log p 

log (-logp) 



Note that, contrary to the fiavour-less case, the dilaton goes to — oo in the IR and not to a 
constant. Moreover, the Hamiltonian (B.9) for Nf 7^ is not simply related to the dilaton. 



(B.52) 



C An analytic solution for Nf = 2Nc 

Here we summarise some results about an analytic solution that is obtained as the rotation 
from the one with Nf = 2Nc. We will use as "seed" solution the solution quoted in eqs.(4.7)- 
(4.8) of the paper [H]. Notice that these solutions are qualitatively different from those studied 
in the main text, for example, the fibrations in the metric are absent, and the function Q{p) 
is constant. The functions P{p),Q{p) are known exactly in this case and given by 



P 



9Nr 



+ ce3 



Q 



3Nr 



cosh r = 1. 



(C.53) 



We will keep only the upper sign in the following. In this case the radial coordinate ranges 
in the whole real axis. For large values of p the geometry asymptotes to the conifold, while 
in the IR, for p — t- — 00, the geometry is the one for Nf = 2Nc with a linear dilaton, obtained 



in eqs.(4.22)-(4.23) of the paper [7] (for the value of 3^ = 4). Using (3.19), we find functions 



47 



of the background read 
,2, _ 



4p/3 



4p/3 



+ -e 
4 6 



4p/3 



-40^[1 + — ][1 + — ^ ^ — £ 1 



2c 



2c 



h 



-20 



2c2 
a = b 



0. 



(C.54) 



The warp factor h for large values of the radial coordinate (p — oo) has the following 
expansion: 



h 



'20c 



27NI 
128c3' 



'16p/3 _^ (^(^g-20p/3^) 



(C.55) 



Asymptotically in p — )■ oo, we have that the dilaton is constant and the metric and fluxes 
read 



ds 
B3 



10 



,2p/3+</.„ 



'-dp^ + -(d(^2 + dfi2) + -(dz/' + cosMv? + cosMy^)^ 



4 

6-^3, 



3Nc ' 3 ' ' 4:'^'''^ ' ' 6' 

sin 9d6 Ad(p + sin 6d6 Ad(f A {d^p + cos 6d(p + cos 9dip) . 

i^5^(l + *io)^^e^^/Sol4Adp. 



(C.56) 



On the other hand, when we explore the geometry near the IR (p — ?■ —00 where there is 
a singularity), we will find that the whole background resembles that in eqs.(4.22)-(4.23) 
of the paper [7] for the value of 3^ = 4. The reason for this is that when going to the far 



IR, the rotation is "undone", according to eq.(2.10) and the dilaton vanishing rapidly as 



I p-->— 00 



Notice that contrary to the Nf = 2Nc solution discussed in the main text, this solution 
preserves the f/(l)_R symmetry, while it does not preserve any Z2 symmetry discussed in 
section |5| Further studies of these special solutions, and their field theory interpretations, 
are left for future work. 



D Flavoured fivebrane solution in the variables of [5j 

For convenience of the reader, in this appendix we write the fivebrane solution in the variables 
used in In particular, the complete solution may be written in terms of two function^ 
c, /, obeying two coupled first order differential equations. The variable t used here is related 

^^The function c in [S] is not to be confused with the parameter c used elsewhere in this paper. 
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to that used in the main text as t = 2p. Following [5] we write the solution for NS fivebranes, 
with the metric in string frame: 



N, 



d4r = dx^+i + ^ds^ 



4 



(D.57) 



tanhr^? + ^^ + ^? + ^^) + 2^(^^^^ + ^^^^^ 



- 1- 



t 



2Nr Vtanht 



l\{el + el-el-el) 



.2*0 



/ 



1/2 



AT,. 



where 



sinh t 



(es + ^3) A 



+ 1- 



c + 



2N, 



ei A 62 + 1 



Nf \ (tcotht - 1) 



ei A 62 + 1 - 



N 



t 



2Nr sinht 



2Nr 



sinht 



dt A (61 A 61 + 62 A 62) 



(D.58) 

(D.59) 

[ei A 62 + 61 A 62) 
(D.60) 



61 = d9i , 62 = — sin6'id0i , A3 = cos^id0i , 

6i+i62 = 6-'^(d^2 + isin^2d(/>2) , 63 = d^/; + COS ^2 • (D.61) 

The SU (2) left-invariant one-forms 6j obey d6i = —62 A 63 and cyclic permutations. The 



functions c{t) and f{t) appearing in (D.57) obey the equations 



/' = 4 sinh^ t c 

c' = — [c^ sinh^ t — {t cosh t — sinh t)^ 



N 



f 



(D.62) 
(D.63) 



f ' 2N, 

where the primes denote derivatives with respect to t. This system is equivalent to the 
second order equation 

2 



4//" - /'(/' + 8/ coth t) + 16 sinh^ t 



1 - 



— f (sinht-tcosht)2 + — f/ 
2NJ ^ ' 2Nr 



.(D.64) 



Setting Nf = 0, the ansatz and the differential equations reduce exactly to those in [5]. 
Notice the three-form can be written as 



Nf 

— ^^(^3 + ^3) A 61 A 62 + closed 



In these variables it is simple to study the Z2 action generated by [33] 
X: 91^02, 01^02, Hs->-Hs, Fs^-Fs, 



(D.65) 



(D.66) 
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where the action on the three-forms is generated by the center of the SL{2, Z) action. When 
Nf = 0, the three-form H3 is invariant under X, while the metric is not invariant. Therefore 
the full solution breaks this Z2 symmetry |1J. On the other hand, when Nf is non-zero, 



is not invariant under X. In particular, the non-closed part of in (D.65) is not invariant, 
while the closed part is still invariant. In the particular case that Nf = 2Nc the metric 
becomes manifestly invariant under X, however the non-closed part of if 3 is not invariant. 
In this case the three-form becomes simply 

TV 

H3 = ^(es + Aa) A[eiAe2 + eiAe2]. (D.67) 

This is manifestly not invariant under X, however it is clearly invariant under the swap of 
two-spheres, without change of sign of H^. Therefore, for Nf = 2Nc the full solution is 
invariant under ^ ^2 , 0i ^ 02 • 

E More on the flavoured quiver 

In this appendix we make some comments on the quiver gauge theory discussed in section 



4.2.1 of the main body of the paper, see eq.(4.41). In particular, here we will study Seiberg 



dualities and beta functions. The quiver may be written as 

SU{N,) X SU{N2) X 5f/(iV^/2)fl„ (E.68) 
and is described in Figure |6] in the main text. 

E.l Going to the IR 

We start performing Seiberg dualities. After z Seiberg dualities, the quiver will look like 

SU [fii^lffi + + (1 - z)iVi] X SU [^(i±lHf _ + (1 + z)N2\ X ^f/(iV;/2)fl„ 

(E.69) 

where the flavour group is untouched. Assuming that the anomalous dimensions of the Ai, Bi 
fields are 'Ja,b = ~2' (^i^h no assumptions made on the dimensions of the quark superfields), 
the beta functions of each group, at each step of the cascade are 

/3i = {-ir3[N, -N,~ z^] - ^(1 - 7,) 

Nf Nf 

/32 = {-ir'3[m -N,-z^]-^{l- 7,) (E.70) 
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and the change of each group in a given step of the cascade i 

AN, = ' [2N, - 2N2 - -f{2z - 1)] 

AN2 = —^^[2N, - 2N2 - -^{2z - 1)]. (E.71) 

Notice that for odd steps, the group 2 does not change while for even one the group 1 does 
not change. If we use as a criteria for the existence of a cascade down the flow, the fact that 
the changes in groups must be positive, whenever we find a negative AiVj means that we 
must stop the duality cascade. This bounds the number of dualities performed, by a critical 
number given by 

^ ...[i . ^(^1 (E.72) 

where with Int we mean the integer part. Following the discussion in the main part of the 
paper, for simplicity let us consider the case when 

N, = [k + l]N,, N2 = kN,, Nf = ^, z, = lnt\^^]. (E.73) 

A I Z - 

Notice that for < A < 00, we can have Nf small or large compared to Ni,N2. Also, notice 
that A must be a rational number. For example, if we restrict the attention in the interval 
Nc < Nf < 2Nc, we have 2 < A < 4. This in turn bounds the maximum number of Seiberg 
dualities to z^: < 3. It may be instructive to consider the cases A = 2,3,4, to notice that 
the pattern is that after one or two dualities, the quiver comes back to itself, or is in a 
distribution of colours and flavours that does not permit further Seiberg dualities. Although 
from the field theory point of view, we could consider the situation in which A ^ 1, in which 
we can do many Seiberg dualities, our string theory backgrounds certainly do not require 
this. In fact, we typically have A ~ Notice that, on the contrary, in the situation 

studied in [21] we have that ^^J,^ — ?■ 0. This gives the possibility of performing many 



Ari-Ar2 

Seiberg dualities, and in fact z^ diverges there. 
E.2 Going to the UV 

We will quote here similar formulas to the ones above, for the flow towards the UV. As 
usual, the idea is to Seiberg dualise the strongly coupled group. The first point is that after 
s Seiberg dualities (we distinguish with s the dualities to the UV while with z those to the 



^^The change of the group is defined as Rank[before]-Rank[after], and is typically positive indicating a 
decrease of the group and degrees of freedom under Seiberg dualities down the cascade. 
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IR), we have a quiver 
SU 



1)^/ + sN, + (1 - s)N2] X SU r^ii±ll^ _ ^ ^1 ^ 1 ^ sUiNf/2)a„ 



4 

(E.74) 

The beta functions at each step are 

/3i = (-ir^3[iV2-iVi-.^]-^(l-7,), (E.75) 

/32 = (-ir3[iV2-iVi-.^]-^(l-7,). (E.76) 

Notice that the beta function for the sum of the gauge couphngs is negative (assuming 
7q < 1), independently of the number of Seiberg duahties. This imphes that at some point 
in the UV at least one of the couplings will diverge, that is, there is a Landau pole. This 
behaviour was indeed observed in the solution of [HI], which was proposed to be dual to the 
quiver we are discussing, but at the origin of moduli space, i.e. not in a Higgsed phase. 
The change in groups is 

1 + (-1Y Nf 
AiVi = ^ [2N,-2N,-^i2s-l)], (E.77) 

AiV2 = \ ^ [2N2-2N,-^i2s-l)] (E.78) 

Now, going to the UV, one would in principle keep on doing Seiberg dualities unless ANi > 0, 
which would mean that there is a decrease in the number of degrees of freedom. This implies 
that 

ANi > ^ s < - - (E.79) 

2 Nf ^ ^ 



Coming back to eq.(E.73), this implies that for s < we have to stop dualising. For 
any s > 0, we have that the inequality below cannot be satisfied if M > 0, Nf > 0, s > 1, 
so, if we can do one Seiberg duality, we can do as many as we want. As we can see, this 
quiver has non-stopping Seiberg duahties to the UV. In the case of [31], one can perform 
many Seiberg dualities before reaching a "duality wall". This occurs at some finite energy 
scale, where the Seiberg duahties "accumulate" , and the number of degrees of freedom (for 
example measured by the central charge) diverge. 



F The holographic central charge 

In this appendix we compute the "holographic central charge" for our backgrounds. To 
define this quantity, the idea is to reduce the ten-dimensional system to a five-dimensional 
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sigma model coupled to scalars, following for example the procedure described in [36] . Then 
one uses results derived in |37], to show that such quantity is monotonic, and stationary 
at AdS points. Thus it may be interpreted as a central function, measuring the degrees of 
freedom of the field theory. Defining the functions 

16 



and reducing our background in eq.(2.11) to a five- dimensional gravity theory coupled to a 
sigma model of scalars, we find (in the Einstein frame of the five-dimensional theory) that 
the metric reads 

d4^ = i^i/3[d43 + /3d/]. (F.81) 

One can define a quantity, that may be identify with the central charge (more precisely, 
central function) as 

c oc (F.82) 
Specifying this to our background we find 



C oc /i2e60+4fc+2/.+2g 



i + (V + 4/i' + 2A;' + 60')! ^- (F.83) 
h -I 



We will evaluate this expression at large radius. Using the UV expansions of the various 
functions, we see the term in square brackets is a numerical coefficient, which we will ignore. 
The factors of c and of Qs cancel between the exponentials and the warp factor, and we pick 



up precisely the coefficient in the numerator of the warp factor in (4.35), namely we get 



c {nf + nf for r oo. (F.84) 

This is the expected behaviour for the theory SU {N^ + n + nf) x SU {Nf/2 + n + nf), in 
the limit that n + nj ^ and n + n/ ^ Nf. Of course the number of D3 branes is running, 
and using this, we have 

c ~ (z/r^ + (iVc - Nf/2f log rf for r ^ oo. (F.85) 

The result for the unfiavoured Klebanov-Strassler theory [v = 0, Nf = 0) was obtained in 
[SB] . We see that in our case the dominating degrees of freedom in the UV are the n f source 
D3 branes. The logarithmic growth is related to the cascading behaviour, whereas the much 
more rapid power-law growth that we see is presumably due to the Higgsing. 

Note that H' does not vanish anywhere, and correspondingly the central function c is 
monotonically increasing up to the UV. In particular, our solution does not display a duality 
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wall, as was observed in the solution of [39], where a divergence of c at a finite value of the 
radial coordinate was shown to exist. If there is a Landau pole in the field theory, it will 
occur at infinite energy, where the number of degrees of freedom measured by the quantity 
c defined above, also diverges. 
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